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Abstract 

We solve exactly the general one- dimensional 0(A)-invariant spin model 
taking values in the sphere 5^"^, with nearest-neighbor interactions, in finite 
volume with periodic boundary conditions, by an expansion in hyperspherical 
harmonics. The possible continuum limits are discussed for a general one- 
parameter family of interactions, and an infinite number of universality classes 
is found. For these classes we compute the finite-size-scaling functions and 
the leading corrections to finite-size scaling. A special two-parameter family 
of interactions (which includes the mixed isovector/isotensor model) is also 
treated, and no additional universality classes appear. In the appendices we 
give new formulae for the Clebsch-Gordan coefficients and 6-j symbols of the 
0{N) group, and some new generalizations of the Poisson summation formula; 
these may be of independent interest. 

KEY WORDS: One-dimensional, cr-model, iV-vector model, RP^~^ model, mixed 
isovector/isotensor model, continuum limit, universality classes, finite-size scaling, 
hyperspherical harmonics. 
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1 Introduction 



The purpose of this paper is to study the continuum hmits and finite-size-scahng 
functions in a general class of one-dimensional 0(iV)-invariant spin models (also 
called nonlinear cr-models). Despite the relatively trivial nature of physics in one 
dimension, this exercise is interesting for several reasons: 

1) Two-dimensional nonlinear cr-models are of direct interest in condensed- 
matter physics, and they are of indirect interest in elementary-particle physics be- 
cause they share with four-dimensional gauge theories the property of perturbative 
asymptotic freedom [1,2,3,4]. In particular, recent work [5,6,7], combining Monte 
Carlo simulations and heuristic analytic arguments, has given evidence for the ex- 
istence of new universality classes for the two-dimensional 0(iV)-invariant lattice 
cr-model with mixed isovector/isotensor action. The present work was motivated 
by the idea of investigating the occurrence of analogous universality classes in the 
one-dimensional case, where an exact analytic treatment is possible.^ 

2) A second motivation was to perform the computation of an exact finite- 
size-scaling function (as well as the leading correction to it) for a non-trivial spin 
model. Finite-size scaling has become increasingly important in the analysis of 
Monte Carlo data [8,9]. (For example, the functions derived in this paper can be 
used for comparison in the multigrid Monte Carlo study of the one-dimensional 
0(4)-symmetric non-linear cr-model [10].) Moreover, finite-size scaling is the basis 
of an important new method for extrapolation of finite- volume Monte Carlo data to 
infinite volume [11,12,13]. It is also useful to know something about the corrections 
to finite-size scaling. In particular, in the new methods for extrapolation to infinite 
volume, it is crucial to understand the corrections to finite-size scaling because they 
induce systematic errors in the extrapolation. 

3) Finally, our solution method makes use of the functions defined by the gener- 
alization of the usual spherical harmonics to the iV-dimensional unit sphere S^~^, 
which we call hyperspherical harmonics. Although these functions are well known 
[14,15,16,17,18,19], we were unable to find any convenient list of their properties 
in the literature, and therefore we thought that it would be useful to make a com- 
pendium of the relevant properties and formulae. In particular, we were unable 
to find the Clebsch-Gordan coefficients anywhere in the literature (although they 
too are probably known). Using the representation of hyperspherical harmonics 
as completely symmetric and traceless tensors^, the computation of the Clebsch- 
Gordan coefficients is a straightforward combinatoric exercise. Indeed, we can go 
further and compute many of the 6-j symbols. We believe that hyperspherical har- 
monics constitute the most efficient approach to the derivation of high-temperature 
expansions for 0(iV)-invariant spin models taking values in S^~^. Indeed, they 

^We thank Erhard Seller for the suggestion to do this. 

^This representation is of course well known, but it is not (so far as we know) employed in any 
of the standard treatises on hyperspherical harmonics. As we shall show here, this representation is 
an extremely convenient one; one of the purposes of this paper is to make some advertising on its 
behalf. 
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have been used for this purpose by the King's CoUege group [20,21,22,23,24] and 
others [25,26,27]; but the methods were cumbersome, in part due to the lack of 
convenient expressions for the Clebsch-Gordan coefficients. In addition to the work 
reported here, we are now using these methods to extend various high-temperature 
expansions for two- and i/iree-dimensional 0(iV)-invariant and ?7(iV)-invariant spin 
models [28]. 

This paper is organized as follows: The hyperspherical harmonics are intro- 
duced in Section 2, where we also explain how they are used in the expansion of the 
Gibbs weight exp(— In Section 3 we give the exact solution for the general one- 
dimensional S^~^ cr-model in finite volume, as well as its infinite- volume limit. All 
expressions are written in terms of the normalized expansion coefficients ujv,/ (which 
generalize the well-known v = tanh J for the Ising case N = 1). In Section 4.1 we 
discuss in detail the possible continuum limits for one-parameter Hamiltonians by 
performing the large- J (i.e. low-temperature) expansion of ujv,/( J), and we show the 
appearance of infinite families of universality classes. The finite-size-scaling func- 
tions and the corresponding corrections to finite-size scaling are given in Sections 
4.2 and 4.3, respectively. Finally, in Section 5, we analyze a class of two-parameter 
Hamiltonians — which includes, among others, the mixed isovector/isotensor model 
studied in [5,6,7] — and we show that no additional universality classes appear be- 
yond the ones already found in Section 4.1. In Appendix A we provide proofs of 
various properties of the hyperspherical harmonics, including the Clebsch-Gordan 
coefficients and some of the 6-j symbols. In Appendix B we analyze the finite-size- 
scaling functions for a one-parameter family of universality classes that includes 
those of the mixed isovector/isotensor model; this analysis is based on generalized 
Poisson summation formulae applied to some generalized theta functions. (We think 
these formulae may be of independent interest; as far as we know they are new.) In 
Appendix C we study the limit iV — )• oo of the finite-size-scaling functions for the 
standard iV-vector universality class. 

2 Hyperspherical Harmonics 

The purpose of this section is to introduce the hyperspherical harmonics that 
will give the basis for expanding the Gibbs weight for our spin models. From 
the mathematical point of view this is connected with doing harmonic analysis on 
the unit sphere S^~^ C R acted on transitively by the compact connected Lie 
group SO(N) [30]. More precisely, let us consider: 



<N-1 



• Re SO{N) 



• the normalized rotation-invariant measure c/0(<t) on S 



<N-1 



• the (complex- valued) square-integrable functions f ^ 
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• the unitary representation T{K) of SO{N) on £"^{3^ ^) defined by (T(i?)/)(<T) 

Then, we want to find an orthogonal Hilbert-space decomposition of L'^{S^~^) into 
subspaces such that the representation T(R) restricted to each subspace is irre- 
ducible. The needed decomposition turns out to be precisely the decomposition of 
L'^{S^~^) into eigenspaces of the Laplace- Beltrami operator jC = £5^-1."^''* In fact, 
it can be proved^ that: 

(a) The eigenvalues^ of C are 

Ajv^fc = A;(iV + A;-2) > 0, (2.1) 
where A; = 0,l,2,.... The corresponding eigenspace Ej^^k has dimension'' 

^ N + 2k-2 T(N + k-2) , , 

M^,, ^ dim E^,, = V^^^TTT ^^-^^ 

and can be given several equivalent descriptions: 



^The Laplace-Beltrami operator on ^ can be defined as follows: Define on the vector 
fields ( "angular-momentum operators") 

r«/3 - ( r"— tP — 

Then the restriction to 5'^""'^ of each L"'' is a vector field on 5'^""'^, and 

l<a<p<N 

is the Laplace-Beltrami operator on S^~^ . 

■*We remark that, for N >?>, C generates the algebra D(5'^~"'^) of 5'0(#)-invariant differential 
operators on 5'^""'^. (For N = 2 this is not the case, because d/dO is an 5'0(2)-invariant differential 
operator not belonging to the algebra generated by C But if we consider differential operators 
invariant under 0(N) instead of SO(N), then the assertion is true also for N = 2.) 

^See [30], Theorem 3.1 (pp. 17-19). 

®Note that our C is the negative of the usual Laplacian, i.e. it is a pos«<«i)e-semidefinite operator. 

^For a proof see [30], Exercise A.5(i) (pp. 74, 552) and [15], Lemma 3 (p. 4). See also Appendix 
A.l below. Usually we are interested in the case > 3, for which formula (2.2) is unambiguous. 
But (2.2) is also valid for # = 1,2, if it is interpreted as an analytic (in fact polynomial) function 
of N for each fixed integer k > 0. Thus, for N = I and # = 2 we have 

A-,, . = ,„ £±|i^ M±izi) = [ 1 to, . = 0.1 

N^i k\ T{N - 1) [0 for A; > 2 

and 

N + 2k-2 T(N + k-2) r 1 for = 



M2,k = dim£'2,fc = lim r^/Ar i\ - 1 o ^ ; \ 1 

N^2 k\ T(N — 1) [2 for fc > 1 

Note also that Wat = 1 and Wat 1 = # for all N. 
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(i) -Ejv,fc consists of the restrictions to S^~^ of the harmonic polynomials of 
degree k on (namely, the homogeneous polynomials of degree k that 
satisfy Laplace's equation on R ). 

(ii) En k is spanned by the functions /(<t) = (a • cr)'^ with a G and 
Ei, a] = 0. 

(iii) -Ejv,fc is spanned by the completely symmetric and traceless tensors 
y)^\"'"*(<T) of rank k, as the indices cti, a2, • • • , ctfc range over the 
allowable values.® (These tensors are described in more detail below.) 
Of course, since in general > dim£^jv,fc, the Y^\""''{(t) form an over- 
complete set. 

(b) Each eigenspace -Ejv,fc is left invariant by T(R). Moreover, for iV > 3 the 
representation T(R) t-Ejv,fc of SO(N) is irreducible.^ 

oo 

(c) L'^{S^~^) = -Ejv,fc (orthogonal Hilbert space decomposition). 

k=0 

To make all this concrete, we can write: 



I^N,k ((t"i • • • a"* - Traces) 



(2.3) 



where a G ^, "Traces" is such that y)^\"'"*(<T) is completely symmetric and 
traceless (namely"^° (5„^„^y)^^^"'"*(<T) = for any i ^ j), and 



Explicit examples are: 



1 



klTif 



1/2 



(2.4) 



(2.5) 



(2.6) 



N 



(2.7) 



/iV(iV + 2)(iV + 4) 
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®The functions (a • a)^ used in description (ii) above are linear combinations of the Y's (see 
eq. (2.3)), namely [a- erf = a„j . . .a„^ ^jv.V ('^)- ^he condition Ya=i '^i = ^ ensures 

that the "Traces" in (2.3) make no contribution. 

^For N = 2 the group is abelian, and the spin-A; representation for k > I decomposes into the 
two irreducible representations e^'*^^. However, if we consider 0{N) rather than SO{N), then the 
representation is irreducible also for N = 2. 

^"^The usual summation convention will be used in this paper from now on. 
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„ct „fJ „-) 

a a a 



1 



iV + 2 



(2.8) 



/iV(iV + 2)(iV + 4)(iV + 6) 



24 



+ 



1 



iV + 4 



1 



(iV + 2)(iV + 4) 



6"^^a"'a^ + 5 permutations 



^al3 pS _|_ pi^pS _|_ pS ^I3f 



(2.9) 



[The general formula is given in equation (A. 17).] We note that for iV = 3 the 
y's are linear combinations of the usual spherical harmonics, and (MvaE^^k = 2k -\- 
1. Similarly, for N = 2 the Y^s are linear combinations of cos k6 and sinA;^ (or 
equivalently of e^'*^^), and dim£^2,fc = 2 for A; > 1. For iV = 1, Yiv,/; vanishes for 
A; > 2, while Ki^o = 1 and Yi^i = a. 

The normalization fJ,N,k is chosen so that the following orthogonality relation 
holds (see Appendix A. 2): 



(2.10) 



where j^i^--"*''^!---'^* jg ^]-^g unique orthogonal projector onto the space of completely 



symmetric and traceless tensors of rank A;, defined by the following properties (see 
Appendix A. 3): 

1. complete symmetry in the indices a, and in the indices fi 

2. symmetry under the total exchange fii for all i 



rpai...af^ 



for any completely symmetric and traceless 



tensor Tjv,fc 
As special cases of condition 4 we have 



LN,k 



and 



N,k ^ N,k i^) 



J- AT I- J- I 

For example we have: 

pf3 
1 



T«i«2;/?i/?2 



'JV,3 



N 



pict2 ^Pifh 



(2.11) 
(2.12) 

(2.13) 
(2.14) 



^ai/?l^a2/?2^a3/?3 _j_ ^"ife ^a2/?3 ^"s/?! _j_ ^ai/?3 ^a2/?l ^"sfe 
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_j_ ^ctiPi ^a2fh ^cto,f32 _j_ ^ctiPz ^a2P2 ^ctzPi _|_ ^a2/?i ^"s/^s 

-'- ^^0.10.2 f ^PlP2 ^CtsPs _j_ ^/3l/33^a3/?2 _j_ ^/?2/?3^a3/?l 



3(iV + 2) 

_j_ ^aia3 ^^/?l/?2^a2/?3 _j_ ^PlPs ^0:2^2 _j_ ^/?2/?3 ^a2/?l ' 
_j_ ^a2a3 ^^/?l/?2^ai/?3 _j_ ^/?l/?3 ^"ife _j_ ^/?2/?3 ^"ift ' 

(2.15) 

[The general formula is given in equation (A. 27).] The trace of this operator is given 
by [see (A.35)/(A.36)] 

= ^N,k = dim E^,, , (2.16) 

as of course it must be. We remark that Y^^kif^) ■ YN,k{f^) = ^jvfc ""*(''')-^JV 
is independent of cr [by 0(N) invariance], and hence 

YN,ki(T) ■ YN,kif^) = ^N,k (2.17) 

by (2.10) and (2.16). 

As stated in the theorem given at the beginning of this section, the hyperspheri- 

cal harmonics are a complete set of functions on L'^{S^~^). Thus any function /(<t) 
can be expanded as 

oo 

fi^r) = E fk'-"' Y^^-^ia) (2.18) 

k=0 

where 

fk'-"" = I d^l{T) fir) Y-]--{t) . (2.19) 

For smooth functions this expansion converges very fast. Indeed, if /(<t) is infinitely 
differentiable, then, for A; — )• oo, the coefficients of the expansion go to zero faster 
than any inverse power of k (see Appendix A. 4). 

The completeness of the hyperspherical harmonics can be expressed through the 
relation^^ 

oo 

Y: Y^:r''{a)Y^:r''{T) = S{a,T) (2.20) 

k=0 

where the (5-function is defined with respect to the measure c/0(<t). 

Finally, let us consider an invariant function of two "spins" <t, r G S^~^, i.e. a 
function of <t-t.^^ We want now to compute its expansion in terms of hyperspherical 
harmonics. Using Schur's lemma (see Appendix A. 4) we can write 

oo 

f(a -t) = FN,k liv,fc(o-) • FAr,fc(r) . (2.21) 
fc=o 



^^Note that the normalization here follows directly from the one defined for (2.10). 

^^For N = 2 there are functions of cr, r which are 5'0(2)-invariant [but not 0(2)-invariant] and 
are not functions of <t • r: namely, they can depend also on cr x t. We are not interested in such 
functions. 



8 



We can drop the "Traces" terms of either one of the Y^s in the scalar product 
above, since the other Y is traceless. Also, since the scalar product is rotationally 
invariant, we can rotate <t to w = (1, 0, . . . , 0) and correspondingly rotate r to 
some p with cr ■ t = w ■ p = . In this way we obtain 

Fjv,fc(cr) • Fjv,fc(r) = Fjv,fc(w) • Fjv,fc(p) 

= fiN,kw"^...w"^Y^^f"'ip) = fiN,kY^fip). (2.22) 

Now Y]^'i^(p) can be expressed in terms of Gegenbauer polynomials^^ (this corre- 
sponds to the relation between Yiq and Legendre polynomials for the usual spherical 
harmonics) as^"* (see Appendix A. 2) 

and therefore 

In particular, for w = (1, 0, . . . , 0), we have 

F-,Hw) = ^ . (2.25) 

From equation (2.21), using the orthogonality relations, the rotational invariance of 
the measure, equation (2.17) and (2.24), we get 



FN,k = I dn(p) f{/) / . (2.26) 



Now the integrand depends only on and we can integrate out the other coordi- 
nates. We finally get 

Fn, = %i / dt{i- er-^y^ m %^ , (2.27) 

J-l C\' (1) 



i^See [31], pp. 1029-1031. 

"'^■^For N = 2 this relation is singular, since C^(x) = 0. This singularity is due simply to the 
normalization convention of the Gegenbauer polynomials, and indeed the limit # ^ 2 is well- 
defined. The result is simply 

Cf/'-^(cos^) Tfc(cos^) 
lim " " = cos k6 = , , . 

cf''-^(i) nil) 

where Tk(0) are the Chebyshev polynomials of the first kind (see [31], formulae 8.934.4 (p. 1030) 
and 8.940.1 (p. 1032)). 
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where Sn is the surface area of the iV-dimensional unit sphere: 

From the general properties of the hyperspherical harmonics we can derive the 
foUowing properties of the coefficients -Fjv,fc (for the proofs of properties 1 and 2, see 
Appendix A. 4): 

1. If f(t) is positive^^ for t G [-1, 1], then |-Fjv,fc| < -Fjv,o for aU k ^ 0. 

2. If f(t) is smooth (i.e. C°°), then hm k^F^^k = for every n. 

3. If f{t) = t', then the integral in (2.27) can be performed explicitly^^ and the 
coefficients -Fjv,fc are given by 

f r(f)r(/ + i) 

/;\ 7 r 7 r if A; + / is even and k < I 

F^], = 2'r(^) r(i^ + l) - (2.29) 

V otherwise 

and are, in particular, always nonnegative. It immediately follows that for a 
generic function 

oo 

m = E /'^'' (2-30) 
the coefficients -Fjv,fc are given by 



= J2fi F^]k ■ (2.31) 



.(0 
i=k 

Therefore, if all the coefficients /; are nonnegative, then so are the -Fjv,fc. 



In particular, using (2.27) or (2.31) it is possible to compute the coefficients -Fjv,fc 
for the functions exp [J (<t • r)] and exp (<t • t)^| . In the ffist case we obtain 



1— — 

FN,k = r(|) (I) ' /^+,_i(J) (2.32) 
where is the modified Bessel function}'^ ; in the second case the integration gives 



r(f) r(^) fJ^^/^ ^(k + i^Nj\ „ 

^ ^ 1-^1 ! ^ H ! ~ for even k 



^-'^ = 1 Ar(f + ^) V2; 2 "^-2'2y (2.33) 

I for odd k 



"'^^More precisely, it suffices that / be nonnegative and not almost- everywhere-vamshmg . 
i^See [31], formula 7.311.2, p. 826. 

^^In particular, for N = 1 (the Ising model) we get Fi o = cosh J and Fi i = sinh J, and there- 
fore the formulae in the following sections will be written in terms of the usual high-temperature 



expansion parameter vi i = Fi i/Fi o = tanh J. 
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where iFi is the confluent (degenerate) hypergeometric function}^ These two ex- 
pansions wiU be used in the next section. 

Let us now compute the Clebsch-Gordan coefficients. In general we can write 

^JV,fc K'^) ^N,l [f^) - 1^ ^N; k,l,m ^N,m [f^ ) [2.64:) 

ra 

Using the orthogonality relations (2.10) we obtain 

C/Mi'"-'""-""" = /rfO(^) (2.35) 
This integral can be computed explicitly. We get (see Appendix A. 5) 

-,ai...au]fh-Pi]li-lm _ l^'N,k I^N,l I^N,m k\ l\ ml 

ik+j)\ Zljlhl 

jai...a,,]ai...aihi...hh ...ft ;6i ...fefeCi ...Cj ^71 ...7™ ;ci ...Cj ai ...a, z^, „„x 
^ ^N,k ^NJ ^N,m [^,.00) 



*^JV; k,Lm 



if |/ — A;| < m < / + A; and A;-|-/-|-m is even, with z = (m + A; — /)/2, j = (m + Z — A;)/2, 
h = (Z+A; — m)/2, and vanishes otherwise. (Of course we are considering A;, /, m > 0.) 
In the following we will be interested in the scalar quantity 

^N;k,l,m ~ ^N; k,l,m ' C]\f- k,l,m ■ (2.37) 

The general formula is reported in Appendix A. 5 [see (A. 63)]. A particularly simple 
case is m = / + A;: 

2 2 

^2 xr ^^N,k f^N,l /r, oo\ 

'-JV; k,l,l+k - ^^N,l+k 2 ' (2.38j 

which can be obtained directly from (2.36), using the properties of the /jv,fc tensor 
and (2.16). If A; = 1 this gives 



Cl,r,,,r = ^{^^1 ^1- (2-39) 



It follows immediately from (2.35)-(2.37) that C^. ^ i ^ is symmetric in the vari- 
ables A;, / and m. This implies, for example, that 

2 2 

'-N;k,l,l-k - ^N;k,l-k,l " ^N,l ^ (2.4Uj 

[from (2.38)]. It also implies that, for k fixed, it suffices to find C^. ^ ; ^ for I < m < 
I -\- k . Thus, the two coefficients needed (for each /) for the case k = 1 are obtained 
from (2.39). For the case A; = 2, which will be used later on, we have from (2.38) 
that 



r2 N{N + 2){N + l-l) N + l-2\ 

^N;2,l,l+2 - 2(iV + 20 [ I ' ' ^ ^ 



^See [31], pp. 1058-1059. 
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and from (A. 63) we obtain^^ 

.2 _ ^ l{N + 2){N-2){N + l-2) 

''''' - (iV + 20(iV + 2/ - 4) • ^^-^^^ 

Using the completeness relation (2.20), formula (2.35) and (2.17) it is easy to 
verify the identity 

oo 

E CN;k,Wra = J^N,lMN,m. (2.43) 
fc=0 

3 Exact Solution for a Generic h 
3.1 Finite Volume 

In this section we want to discuss the most general O (iV)-invariant cr-model tak- 
ing values in S"^"^, with nearest-neighbor interactions, defined on a one-dimensional 
lattice with L sites and periodic boundary conditions. We consider a Hamiltonian 
of the form 

L-l 

'H{{(t]) = -J2 h(a,-a,+,) (3.1) 
with (Tl = (Tq. Interesting special cases are the N -vector model 

h{(T^ ■ ay) = J (T^ - (Ty (3.2) 

and the RP^~^ model 



h{(T^-(Ty) = ^{(T^-(Tyf . (3.3) 



The coefficients -Fjv,fc have already been evaluated for both of these models [see 
(2.32) and (2.33)]. 

We want to evaluate the following quantities: 



• Partition function: 

L-l 



ZN{h-L) = j'^o- J] e''('^-'^-+i) (3.4) 

x=0 



Spin- A; two-point function (A; = 1,2, ...): 



GN,k{x,h;L) = — ^ ( Fjv,fc(cro) • Fjv,fc(cr^) )l (3.5) 

L-l 

GN,kip,h-L) = ^'""^ GN,k{x,h-L) (3.6) 

x=0 



"'^^Formula (2.42) is potentially ambiguous if # -|- 2/ — 4 = 0, which can happen for {N = 2,1=1) 
and {N = 4, / = 0). In fact C'j^^ 211 = N{N — 1) and C'j^^ ^ ^ ^ = 0; these results can be obtained by 
interpreting (2.42) as an analytic function of N for each fixed /. 
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where < x < L, and p is an integer multiple of 27r/L. Note that the 
normalization A/jv,fc [defined in (2.16)] ensures that GN,k (0, h; L) = 1.^° 

• Susceptibility (= two-point function at zero momentum): 

XN,k(h;L) = GN,k(0,h-L) (3.7) 

• Two-point function at the smallest nonzero momentum: 

• Second-moment correlation length: 

(h; L) = 2sin(vr/i:) W > ^N,k ^g^g^ 

I undefined otherwise 
In all of these formulae, we have used the abbreviations 

L-l 

Va = Y[ dn(a^) (3.10) 

(/(M))l ^ JVa /(M) e-»(^^'» (3.11) 

To compute all these quantities, we expand in terms of the hyperspherical 
harmonics Y^^k^ as described in the previous section: 



oo 



exp[h(a^ ■ ay)] = ^ FN,k(h) Fjv,fc(cr^) • Fjv,fc(crj^) . (3.12) 

k=0 

The integration over T>(t is then immediate using the orthogonality relations (2.10) 
and the integral (2.35). In this way (using also the symmetry of C^. ^ ; ^ in the 
indices / and m) we obtain: 



ZNih;L) = FN,oihf J2 ^N,iVNAhf (3.13) 



^"^Since we are using periodic boundary conditions, G]\[,k {x, h; L) = G]\[,k {L — x, h; L) and there- 
fore GN,k (p, h; L) is real. Indeed, Gat ^ {p, h; L) > for all p because equation (3.6) can be written 
as 

L-l 



GN,k (p, h;L) = ^ 



J2 e'''">iv,fc(^.) 



x = 



by using translational invariance. On the other hand, GN,k (x, h; L) may in some cases be negative 
( "antiferromagnetism" ) . 
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GNA.^.h-L) = E Cl.^KWraVNAhTv^.iKf-^ (3.14) 

Gjv,fc (p, /i; L) = ^^'fl^n £ I ^^JV,/(^)^ 



VN,l{h? - 2(cOSp) VN,l(h)vN,m(h) + UjV,m(/i)2 

where < x < L — 1; here we have defined the normahzed expansion coefficients 

VN,k{h) = -— , (3.16) 

-rjv,o [It-) 

which wiU play a central role in the subsequent analysis. Notice that, because of 
the properties of the coefficients -Fjv,fc discussed in Section 2, we have |ujv,fc {h)\ < 1 
for all k ^ 0.^^ Moreover, all these series converge very fast (at least if h is smooth): 
this is because, for k — )• oo, v^^h) goes to zero faster than any power of A;, while 
A/jv,fc ~ k^-^ and [see (2.43)] 

'^N;k,l,m < min(A/'jV,fcA/'jV,/, J^N,kJ^N,m, J^N,mJ^N,l)- (3-17) 

Finally, using (2.43), it is trivial to check that Gjv,fc (0, h; L) = 1 in (3.14). 



3.2 Infinite Volume 

We want now to consider the infinite- volume limit L — )• oo in the expressions 
from Section 3.1, keeping the parameters of h fixed. Since |t^jv,fc(/i)| < 1 for all 
A; 7^ 0, V]y^k(h)^ goes to zero for L — )• oo unless A; = 0. Thus in (3.13), (3.14) 
and (3.15) only the term with / = survives in the infinite- volume limit. Since 

^^The summand in braces in (3.15) is potentially ambiguous in two cases: 

(i) p=0 and VN,i = VN,m; 

(ii) p= IT and VN,i = -VN,m- 

In these cases the correct summand is k i m VN,i{h)^ , as can be seen by going back to (3.14) 
and performing the sum over x. [The same result can be obtained formally by symmetrizing the 
summand in / and m (using = k,m,i)' i-'^- replacing vf^j by (vf^ j - v'^ „^)/2, and then 

treating t;jv,; and viq^m as independent variables for which one can take the limit viq^m ifw,;-] 

^^For the RP^~^ model, or more generally if h {a^ ■ (Ty) is an even function, all the coefficients 
Fn,i (h) [and the corresponding t;jv,; (h)] with / odd are equal to zero (by symmetry). Therefore, 
in the above formulae, only even values of / and m can appear in the sums (except, of course, in 
GN,k(x, h; L) for x = 0). From this and the properties of the quantities k i i-^- that they are 
nonzero only if A; -|- / -|- m is even, it follows that for k odd the spin-A; two-point function vanishes for 
all X ^ Q. Of course, this follows equivalently from the 2'2-gauge-invariance of the model when h is 
an even function. 
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^N-kOm — ^kmJ^N,k, we get the well-knowii results^"^ [29,32] 

ZN{h;L) = FAr,o(/i)^ [l + 0(e-^^)] (3.18) 
GNA^,h;L) = VN,kihP fl + 0(e-^^)l (3.19) 



where 



1 -2{cosp)vN^k[h) + VN,kW 



ji = - min log |ujv fc(/i)| . (3.21) 



In particular we obtain 

1 - VN,k{h) 



1 + VN,k{^) /q oo^ 

XN,k[h]oo) = — (3.22) 



and 

f VN,k{hf'' -f 
Ct'\h;oc) = l-mAh) (3.23) 

[ undefined if ujv,fc < 

Let us notice that in infinite volume the correlation functions are simple expo- 
nentials. In fact, if we define the masses mjv,fc(/i) for A; = 1, 2, . . . by 

( -log VN,k{h) for < ujv,fc < 1 

mN,kih) = \ (3.24) 
[ undefined for — 1 < ujv,fc < 

then, in the usual case^"* in which ujv,fc > 0, the correlation functions are 

GnAx, h;oo) = e-"^.*!"! . (3.25) 

We can also define the exponential correlation length by 

C\^f{h;^) = hm = . (3.26) 

4 A One-Parameter Family of Hamiltonians 

In this section we want to study the continuum limits and finite-size-scaling 
functions in a one-parameter family of interactions of the form 

h{a-T) = Jh{a-T) , (4.1) 



^^The formulae in Section 3.1 are written for x > 0. By translation invariance, we obviously have 
GN,k(x) = GN,k(—x). Therefore, we can obtain formulae valid for all x by systematically replacing 
X by I a; I; we have done that here. 

^■*As we will see in Section 4.1, the case of negative t;jv,fc does not give rise to a valid continuum 
limit. 
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where h is some fixed function. Therefore, -Fjv,fc, ujv,fc and aU the quantities intro- 
duced in the previous sections are now functions of J. As /i is arbitrary it suffices 
to consider the case J > only. Since we are in one dimension, there are no critical 
points at finite J; the only way of obtaining a continuum limit is to take J — )• +00. 
We will do this by obtaining an asymptotic expansion of the coefficients i^jv,fc(<^) 
for large J. Using the general formula (2.27) with f(t) = exp[J/i(t)], the prob- 
lem reduces to expanding the integrand around the absolute maxima of h(t) in the 
interval [—1,1]. 

In Section 4.1 we will study the continuum limit in infinite volume. In Sections 
4.2 and 4.3 we will study the finite-size-scaling limit and the corrections to it. 

The discussion in Section 4.1 of the possible continuum limits will be restricted 
to the case iV > 3, since N = 2 displays different properties (related to the different 
topological structures of the sphere for iV > 3 and iV = 2, and to the fact that 
the only nontrivial normal subgroup of 0(N) for iV > 3 is {i/}, while for N = 2 
there are many others). Although for N = 2 the analysis of possible continuum 
limits is not complete^ it is nevertheless valid for the limits included, and so are the 
finite-size-scaling functions and their corrections.^^ 



4.1 Continuum Limits and Universality Classes for > 3 

4.1.1 Generalities on Continuum Limits 

Consider a sequence ( • )^"^ of infinite- volume lattice models. A continuum limit 
is defined by choosing length rescaling factors S^") — )• 00 and field-strength rescaling 
factors C^), such that the limits^'' 

lim d:lGS;i(s('^)x) (4.2) 
lim d:ls('^)-'GS;i(sH-^p) (4.3) 

exist (in the sense of distributions), where d is the spatial dimension. (For simplicity 
we are considering only the two-point correlation functions.) In other words, a 
continuum distance of x centimeters corresponds to x = S^^^x lattice spacings; and 
conversely, one lattice spacing corresponds to S^"^ centimeters, which tends to 
zero in the limit. 

In our case of a d = 1 nearest-neighbor model, the correlation functions are 
pure exponentials [see (3.19)]; the only parameter is the mass parameter v^]s^\. It 

^®In particular, the discussion following (4.36) does not apply for N = 2. 

^®The case # = 1 is even more trivial, as the only possible function h is h(t) = t. 

^^We use the Fourier-transform convention 



^N,k \^) — 
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is easiest to work in p-space: for any fixed continuum momentum p, the lattice 
momentum p = S^") p tends to zero as n — )• oo , so we can approximate 

cosp ^ 1 - = 1 - S('')"^pV2 • (4.4) 

Thus, the denominator in (3.20) is 

(i-4:i)^ + s('^)"p^4:i. (4.5) 

(Note that S^*^)"^ 0.) 

Now consider a ratio of the correlation function for two different values p, p'. 
If 1 — v^]^\ does not go to zero as n — )• oo at least as fast as S^") ^, then the ratio 

^^nT\p') I ^^nT\p) ^%°T\p) independent of p. This is a physically 

trivial theory (white noise). On the other hand, if 1 — v^^\ goes to zero faster 

than S^") ^, then the limit (if any) will be const/p^, i.e. a massless free field, which 
is ill-defined in dimension d = 1. Therefore, a sensible continuum limit can be 
obtained only when the product (1 — u^^^)S(") tends to a nonzero finite constant 
(which is of course A;-dependent); and this limiting constant is in fact the mass 
"^ivT*^ of the theory. Moreover, in dimension d = 1 li easily follows from (3.19) 

that Cf^\ should likewise tend to a nonzero finite (A;-dependent) constant CwT*^j 
the continuum correlation function is then a massive free field 



— 9 , (cont)'^ 
P + ^N,k 



with mass 

m^-f ^ jim m^^l = jim (1 - 4")) (4.7) 
and field-strength normalization 

^(cont) _ r) ^(cont) (cont) f A Q\ 

Going back to x-space, we have 

G^^T\x) = exp [ - m^-f \x\ ] . (4.9) 

In summary, continuum limits can be obtained from sequences of lattice theories 
in which u^^^ — )• 1 (i.e. m^j^\ — )• 0), and only from such sequences. In particular, 
continuum limits in this sense cannot be obtained from sequences of theories in 
which v^f^\ — )• —1, i.e. antiferromagnetic models with slow decay of correlations. 
As can be seen from (4.9), a continuum limit is uniquely defined by the limiting 
masses and the limiting normalizations CwT*^- Moreover, we shall consider 

two continuum theories which differ only by rescalings of x and the field strengths 
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to be essentially identical. We therefore label the different universality classes by 
the limiting mass ratios, defined as^® 

(cont) (n) 

n^k ^ = lim^. (4.10) 

"''N,k "''N,k 



In the case at hand [Hamiltonians of the family (4.1)], we are considering a 
sequence of theories h = Jh parametrized by J (which plays the role of n). As 
already mentioned, the only possibility for having ujv,fc( J) — > 1 is to let J — )• +oo. 
In the next subsection we will perform an asymptotic expansion of ujv,fc( J) for large 
J, and we will typically find a behavior of the form^^ 

VN,k{J) = l-aN,kA(J) + o(A(J)) , (4.11) 

where the mass scale A(J) and the coefficients ajv,fc will be computed in each case. "^^ 
In this situation, S(J) should clearly be taken to be proportional to A(J)~^, and 
the continuum masses will be 

m^^T^ = aN,k lim A(J) S(J) . (4.12) 

Remark: For some choices of h we will find that the mass parameters ujv,fc(J) 
behave differently according to whether k is even or odd. In such a case we shall 
take A( J) to be of the order of the smallest mass in the theory — which, it turns 
out, is always in the even sector — and we shall write (4.11) only for k even. 
We shall then take S(J) proportional to A(J)~^, and obtain a good continuum 
limit in the even sector. Of course, in the odd sector we have simply white noise 

/ (cont) \ 

4.1.2 Two Simple Cases; iV- Vector and RP^~^ Universality Classes 

Before considering the general case of one-parameter Hamiltonians, let us discuss 
two simple cases of Hamiltonians which generalize, respectively, the iV-vector model 
and the RP^~^ model: 

First simple case: t = +1 is the only absolute maximum of h{t), and h'{l) > 0. 
[This is a subset of what will later be called the Hamiltonians of Type I .] 



^®We choose mjv,2 in the numerator for reasons that will become clear later. 

^®In Section 4.3 we will assume an expansion to the next order [see (4.90)], which will be used to 
compute the corrections to finite-size scaling. In Section 4.1.2 we will explicitly compute such an 
expansion for two simple Hamiltonians (the first case and those belonging to the second case with 
h even): see (4.14) fl". 

Obviously there is some arbitrariness in the definition of A(J): if A(J) is a function satisfying 
C = lim A(J)/A(J) with < C < cxd, then the pair A(J), a^^ k = CaN,k is just as good as the 

pair A(J), aN,k- 
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Starting from (2.27) we first expand the integrand around t = 1 using the relation 



c 



N/2-1 



2Fr[N + k-2,-k- 



N -I l-t 



(4.13) 



where 2-^1(0, c; 2;) is the hypergeometric function.^^ Then, extending the integra- 
tion in t from [—1,1] to [—00, 1], we obtain the asymptotic expansion 



with 

fNiJ) 



FN,k{J) = fN{J) 



l-^ + ^ + 0(J-) 



27rJ/i'(l) 



1/2 



r 



(f ) {'¥ 



l-JV/2 



(4.14) 



(4.15) 



2h'{l) 



1 -1 



(4.16) 



1 j {N + 2k + 1){N + 2k- 1){N + 2k- 3){N + 2k - 5) 



ih'{iy [ 16 

{N + 2k- 1){N + 2k- 3)(iV + 3)(iV + 1) 



(iV + 5)(iV + 3)(iV + l)(iV - 1) 2 (iV + 3)(iV + l)(iV - 1) 
H — r — z s 
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6 



(4.17) 



where Ajv,fc (^0) are the eigenvalues of the Laplace-Beltrami operator on the sphere 
[given in (2.1)], and we have defined 

r = h"{l)/h'{l) (4.18a) 
s = h"'{l)/h'{l) (4.18b) 

For the normalized expansion coefficients ujv,fc( J), we therefore have 



VN,k{J) = 1 



+ 



2h'{l)J Ah'ilfP 



+ 0(J-3) 



where 



(iV + l)r - 1 



(4.19) 

(4.20) 
(4.21) 



^See [31], formulae 9.100 and 9.14.2. 
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The iV-vector model corresponds to h(l) = h (1) = 1, r = s = 0. Notice that in this 
case formulae (4.14)-(4.17) could alternatively have been gotten through a direct 
expansion of the Bessel functions in (2.32). 

Thus, for J — )• +00 all masses [see (3.24)] go to zero as 

mN,k{J) ~ ^N,kHJ) , (4.22) 

where A(J) = l/[2J/i'(l)] is a non-universal scale factor that goes to zero for 
J — )• +00. (Here means that the ratio of the left and right sides tends to 1 as 
J — )• +00.) If we consider the mass ratios defined by"^^ 

K.vJ./)-^^, (4.23) 

we obtain, in the continuum limit, 

7^iv,fc = ^ . (4.24) 

Therefore, all these Hamiltonians give rise to the same continuum limit and belong 
to what we will call the N -vector universality class. 

Second simple case: t = ±1 are the only absolute maxima of h(t) [hence 
h(l) = /i( — 1)], and h' (1) > 0, h' ( — 1) < 0. [This is a subset of what will later be 
called the Hamiltonians of Type II .] 

In this case one must sum the contributions of the two maxima, that is i^jv,fc(<^) = 
Fj^f,(J) + Fj^f,(J). The contribution -Fjyfc coming from t = 1 has already been 

computed. Using the fact that C^^^ = ("1)*^^^^^^ ^(l)? that the 

contribution -Fj^^ coming from t = —1 can be obtained from by replacing 

the derivatives h^'"'^!) with ( — l)"/i(")( — 1) and then multiplying the whole thing by 
( — 1)*^. Thus, keeping only the leading terms, we get 

= 1 - + o(J-') ,4.25, 

^ 2J /i'(i)(i-iv)/2 + |/,/(_i)|(i-Jv)/2 ^ V ; V ; 



for k even, and 



VNk(J) = ^^^^ + 0(J-M (4.26) 



for k odd 
From 

that the odd-spin sector of the theory remains non-critical even at J = -|-oo. [In the 



From these formulae we immediately see that lim ( J)| < 1 for k odd, so 

j^-t-00 



^^We use mjv,2 rather than mjv,i in the numerator in order to facilitate comparison with the 
second simple case below, in which mjv k = +00 for all odd k. 



20 



special case where the function h is evett"^ (as in e.g. the RP^~^ model), we have 
in fact ujv,fc { J) = for k odd, for all J; while for k even we get the same ujv,fc { J) 
as in the first simple case, given by formulae (4.19)-(4.21).] On the other hand, the 
even-spin masses go to zero as 

fnN,k{ J) ~ Ajv,fc A(J) for k even , (4.27) 

where again 

A(J) = ^' (4.28) 

^ ' 2J /i'(l)(l-iV)/2+ |/,/(_l)|(l-iV)/2 ^ ^ 

is a non-universal scale factor that goes to zero for J — )• -|-oo. Thus the limiting 
mass ratio TZ^^k in this case is the same as in the iV-vector universality class for 
even k and is zero for odd k. That is, 

nN,k = ( ^^'2/^^.^ ^ ™ (4.29) 
I for k odd 

All these Hamiltonians belong to the same universality class, which we will call the 
JlpN-i universality class. Notice that the exact Z2 gauge symmetry, which holds 
for the usual RP^~^ Hamiltonian (and more generally whenever h is even), plays 
here no role. Provided that the Hamiltonian has a two-maxima structure at t = ±1 
with h (±1) 7^ 0, the continuum limit will be Z2-gauge-symmetric. For instance, a 
Hamiltonian with h(t) = t'^ -\- a(t — t"^), with |a| < 1 belongs to this universality 
class. 



In summary, we have thus far defined two universality classes: 

(i) the iV-vector universality class, where all the masses go to zero as J — )• -|-oo 
at the same rate and the limiting mass ratio TZ^^k is given by (4.24) for all k; 
and 

(ii) the RP^~^ universality class, where as J — )• -|-oo the even sector displays the 
same behavior as for the iV-vector universality class [i.e. the masses go to zero 
at the same rate with TZ^^k given by (4.24) for all even values of k] while in 
the odd sector the masses either 

(a) do not go to zero [as in the second simple case above] or else 

(b) go to zero at a rate slower than for the even sector [as will occur in some 
examples below], 

and therefore TZ^^k is zero for all odd values of k. 

^^In this case we have FN,kiJ) = for k odd. For k even we have FN,kiJ) = '^^Nki'^)' ^^'-^ 
therefore the coefficients ajv,fc and &jv,fc are given by (4.16) and (4.17), and fN(J) has twice the 
value in (4.15). For r = 1 and s = we obtain the expansion for the RP^~^ model, which can also 
be obtained by direct expansion of the coefficients (2.33). 
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Formulae (4.22) and (4.27) had to be expected on general grounds. Indeed, the 
continuum limit of the iV-vector model (or more generally of any model belong- 
ing to the first simple case above) is simply Brownian motion on S^~^, and the 
generator of Brownian motion is the Laplace-Beltrami operator. "^"^ Thus we expect 
i^N,k{J) ~ A(J)Ajv,fc where A(J) is a non-universal scale factor depending on the 
chosen sequence of lattice Hamiltonians. An analogous discussion applies to the 
JlpN-i cagg; here the continuum limit is Brownian motion on RP^~^ , and thus the 
corresponding masses are related to the eigenvalues of the Laplace-Beltrami opera- 
tor on RP^~^ (which are simply the even-spin eigenvalues of the Laplace-Beltrami 
operator on S^~^). 



4.1.3 General One-Parameter Family 

We want now to address the general problem of studying the limit J — )• -|-oo 
for an arbitrary interaction h; in particular, we want to know whether the two 
universality classes we have just discussed are the only ones which can appear as a 
critical limit of interactions of the form (4.1). As we shall see, the situation is much 
more complicated than this, and in fact an infinite number of universality classes 
appears. 

Let us assume henceforth that h is smooth, and that it has finitely many absolute 
maxima, all of finite order. In particular, suppose it has M absolute maxima on the 
interval [—1, 1] at points ti, . . Am with h(ti) = . . . = h{tM) = h^ax- Let Ui be the 
order of the maximum at t^, i.e. the smallest (nonzero) integer such that h^^'\ti) ^ 0. 
(When ti ^ ±1 the order Ui is of course even and > 2, and h^^'\ti) < 0. When 
= -1 we have h^''^''>(ti) < 0, and when = +1 we have (-iy'^'h^''^''>(ti) < 0.) For 
J — )• -|-oo we have 

M 

F^,k ~ E , (4.30) 

8 = 1 

where Fj^\ is the contribution of the z-th maximum; to leading order in J it is given 
by 



^JV/2-l ^ 

cf/-^(i) 



where 



a; 



r(iV-l)/(2nO if U = ±1 
I 1/n, if U + ±1 



(4.32) 



and Ai is a positive constant, independent of J and A;, given explicitly by 



A. = \ ^ 1 ;:7^ ^ -4. (4.33) 



^■*An arbitrary second-order elliptic differential operator on a manifold M generates a diffusion 
process on M; Brownian motion is the special case in which the generator is the Laplace-Beltrami 
operator. For the general theory of diffusions on a manifold, see e.g. [33, sections 4.1-4.3]. 
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where 

~ _ f 2(^-3)/2 if ^ 

For J — )• +00 the leading contribution comes from those terms with the smaUest 
ctj-; we caU these maxima the principal maxima. Setting a = mincti, we thus have 

i 

F^,, ^ e'^-- A, 'ni^I:: ■ (4-35) 

We want now to know under what conditions the mass mjv,fc(J) tends to zero as 
J — )• +00. For this analysis it is sufficient to use the leading-order expansion (4.35). 
Equivalently we want to see under what conditions ujv,fc(J) — > 1, i.e. when (notice 
that Ca''^~^{t) = 1 and A, > 0) 

E A- = E A- • (4.36) 

Since"^^, for iV > 3, |C^^^ ^(t)| < C^^^ ^(1) when t ^ ±1, this condition cannot be 
satisfied for any k if there is in the sum an i such that ti ^ ±1. Thus the principal 
maxima can only be at 1 or —1. Moreover, if t = —1 appears in the sum, the 
condition can be satisfied only for even k, since C^^^ ^(~1) = ("1)*^^^^^^ ^(l)- 
end up with the following result: 

1. If = 1 is the only principal maximum of /i(t), then 

lim VN,k{J) = 1 (4.37) 

J^ + oo 

for all A; > 1. In this case all correlations become critical. 

2. If ti = ±1 are the only principal maxima of /i(t), then 

/ T\ f 1 for k even , , 

^lim^^^iv,^^) = [,^^^ for A; odd ^^'^S) 

with — 1 < cjv,fc < 1. In this case only the even-spin sector becomes critical. 
In detail, we have 



|/iW(l)|-« + |/iW(-l)|-« ^ ^ 

where n = (iV — l)/2a. 
3. If ti = —1 is the only principal maximum of /i(t), then 

hm v.MJ) = W f'tTf (4.40) 
J^-l-00 L — 1 for K odd 

As in the preceding case, only the even-spin sector becomes critical. 



^See Appendix A. 4. 
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4. If there exists at least one ti ^ ±1 such that = a (i.e. there are principal 
maxima other than ±1), then 

lim VN,k(J) = (4.41) 

with — 1 < < 1 for all k. In this case there is no continuum limit for any 
k. In particular, if there is exactly one principal maximum, and this is a point 
ti ^ ±1, then 



^N^k — • (4.42) 



These results can be understood heuristically: If some ti ^ ±1 contributes at 
leading order to the asymptotic expansion of -Fjv,fc (case 4), then for large J the 
typical configurations have ■ (t^^j^i ti on a significant fraction of the bonds. 
For iV > 3 there are many configurations on each bond with this property (since 
the azimuthal angles are undetermined), and they keep the system disordered even 
at J = +00. In case 1, by contrast, the system orders and thus for J = +oo the 
correlation length becomes infinite. In case 2 the system orders modulo a sign; the 
even-spin correlations are insensitive to the sign and thus display critical behavior, 
while the odd-spin ones remain disordered even at J = -|-oo. In case 3, the system 
develops antiferromagnetic order as J = -|-oo; the even-spin correlations are insen- 
sitive to the sign and thus display critical behavior, while the odd-spin correlations 
have no continuum limit. "^^ In the following we will disregard the theories belonging 
to case 3 [since for the odd-spin sector they do not have a continuum limit, and 
for the even-spin sector they are identical to theories of case 1 with h{t) — )• h{—t)\ 
and to case 4 (since we have proven that they do not exhibit any non-trivial critical 
behavior). The Hamiltonians described in case 1 (respectively case 2) will be called 
Hamiltonians of Type I (respectively Type II). 

To characterize the different universality classes we want now to derive the 
behavior of the masses mjv,fc in the limit J — )• -|-oo. In order to do this, we must carry 
the asymptotic expansion of i^jv,fc(<^) to the first subleading order for the principal 
maxima, and also consider the leading contributions from the non-principal maxima. 
Let us first consider theories of Type I. The relevant expansion for -Fjv,fc is (we 
set = ti = 1) 



IT' ^J^max ) A 7~ 



^®The antiferromagnetic case 3 can be^transformed into the ferromagnetic case 1 by the change of 
variables (fx = {—^Y^'x together with yC{t) = h(—t). The correlation functions then transform as 
GN,k{x, tf; oo) = (— l)*^^GAr,fc(a;, h; oo) and 

/ r- X f GN,k(p,h;(X)) for A; even 

I GN,kip + T^, h; oo) for A; odd 

Thus, case 3 is identical to case 1 for the even-spin correlation functions; and it has no continuum 
limit for the odd-spin correlation functions (since there is no divergence at p = 0). This mapping 
also works in finite volume, provided thai L is even. 
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where 



+ EA- ^Vv. 1^"" ^ (4.43) 



1 r(f±i) 



iv - 1 r 



JV-l' 



2)14 



|/i(-+)(l)| 



(4.44) 



1 - I h^''++^Xl) 

dk+ = \n k + -(N - 3)(N - 1) —- ^ (4.45) 

^.'^ ^ 4^ A ; 2n+ n+ + l M-+)(l) ^ ^ 



leadine term depends now on the relation between a and B = min a,-. We have 



(Note that c+ does not depend on A;, while c//;^+ does.) The first correction to the 
term depends now on the relation be' 

mN,kiJ) ~ AiJ)XN,k + o(J-i/'^+) 

with 

A(J) = c+ J-^/"+ . (4.47) 

(We call the maxima with = /3 the next-to-principal maxima.) Here it should be 
understood that only the dominant term is to be kept: 

(a) If /3 > a + l/n+ , the first term is dominant and the model belongs to the 
iV-vector universality class (4.22). 

(b) If /3 < a + l/n_|_ , then the third term (the term of order 1/ J^~") dominates 
provided that its coefficient is not zero. The coefficient is zero if k is even and 
the only next-to-principal maximum is ti = —1; otherwise the coefficient is 
nonzero. Thus, for k odd the mass is 

^+-^ z:a,=f3 V Cfc (ij/ 

(bl) If the coefficient (for all k) is nonzero, the limiting mass ratio (for any 
k) is 

E_ A- (l 

7^iv,. = 7 ^WttV • (4-49) 



„JV/2- 




qN/2 






"'(to 


„JV/2 
k 


-'(1) 



Clearly there is a multi-parameter family of new universality classes, 
obtainable by varying the {ti,Ai}fl-^ appropriately. 
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(b2) If the coefficient (for k even) is zero, the behavior depends on whether 
the correction o(J~^^~"'^) is larger or smaUer than J~^/"+. If it is smaUer, 
then (for k even) the term A( J) Ajv,fc wiU dominate and therefore, the ra- 
tios 7^jv,fc wiU be those of the RP^~^ universahty class. If the o{ J~^^~"'^ ) 
correction is larger than or equal to J~^l^+ , then a more detailed investi- 
gation is needed. (We note that, also in this last case, the limiting mass 
ratio is zero for k odd just as for the RP^~^ universality class.) 

(c) Finally, if jS = a -\- l/n_|_ , then both terms are of the same order. Again we 
obtain new (multi-parameter) universality classes. In particular, if the only 
next-to-principal maximum is = — 1, we get 

^ _ /Ajv,2/Ajv,fc for A; even fA Kn\ 

^""'^ - \ \N,2l{\N,k + B) for k odd ^^-^^^ 

where < -B < oo is a parameter that interpolates the limiting mass ratio 
between the iV-vector and the RP^~^ universality classes. [Explicitly: B = 
2A/{c+A+).] 

For theories of Type II we set t+ = ti = 1 and t_ = ^2 = —1; the expansion 
of -Fjv,fc is then given by 

« { .4, J- (l - + (-1)' A. J- (l - + 

+ o{J—'^-) + E A. [j-"' + o(J--)] (4.51) 

8 = 3 C/ ^ ( 1 ) I 



where n = n_|_ = n_ = (N — 1) /2a; here dk^+ and c_|_ are given by formulae (4.45) 
and (4.44), and dk^- and c_ can be obtained from the same formulae by simply 
substituting h^'"''^(l) with ( — l)"/i(")( — 1). Defining /3 = min^cti, we obtain that, 

for even values of A;, the masses are given by equation (4.46) with Ai replaced by 
-|- A_ and 

(as before, only the dominant term should be kept) while, for odd values of A;, we 
have 

- 1 < hm ^;^.,,(J) = ^+ ~ < 1 . (4.53) 

Therefore the odd-spin sector of the theory is always non-critical, while the even- 
spin masses go to zero at the same rate for all even k. It follows that the mass ratios 
7^jv,fc = i^N,2/'m-N,k are zero for odd k and nonzero for even k. If/3>a-|-l/n, 
we reobtain the RP^~^ universality class (4.27), while in the other cases an infinite 
number of new universality classes appear. 



See case (b) of the definition of the RP universality class in Section 4.1.2 
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4.1.4 An Example in More Detail 

Finally, to examine more closely the possible universality classes, let us consider 
the special case in which there are only two possible maxima, namely those at 
t = ±1. This generalizes the two cases studied in Section 4.1.2: in that section 
we required h'(±l) ^ [i.e. n± = 1], while now we lift this restriction. Let n_|_ 
and n_ be the orders of the first non- vanishing derivatives at t = +1 and t = —1 
respectively; we will suppose n_|_ > n_ (since we are considering only theories of 
types I and II). Then, from the previous discussion we find four cases: 

For Hamiltonians of Type I, namely n_ < formula (4.46) becomes 



C 



^N,k , [l-(-l)^]A 



^NdJ) ^ -jrr^ + ' jlJ^l - (4.54) 

with a = and fi = a_. 

There are therefore three possibilities: 

(a) If n_ < [{N - l)/(iV + 1)] n+, we get 

mN,k{J) ~ A(^)Ajv,fc (4.55) 
for all A;, with A(J) given by 

A(J) = c+ J-^/"+ . (4.56) 
Therefore, the model belongs to the iV-vector universality class. 

(b) If n+ > n_ > [{N - l)/(iV + 1)] n+, we then have 

rr, ^7^ ~ /^('^)^^'^ for A; even , 

ruNAJ) ~ I 2 (A_M+) J-if'--) for k odd ^^"^^^ 

with A(J) given by (4.56). 

Therefore all masses go to zero as J — )• +oo, but with different rates, so that in 
the limit the odd-spin masses are infinitely larger than the even-spin masses. 
The limiting mass ratios TZ^^k are those of the RP^~^ universality class. 

(c) If n_ = [(N — l)/(iV -|- 1)] n_|_, both terms in (4.54) contribute at the same 
order. We obtain 

^ j HJ) ^N,k for k even 

^^''=~\A(J)(A^.,, + i?) for A; odd ^^'^^^ 

where A(J) is given by (4.56) and 

2 A 

B = (4.59) 
c+ A+ 

is a positive constant. So we get an infinite number of different continuum- 
limit theories, parametrized by B. Notice that < -B < oo; therefore, the 
iV-vector and the RP^~^ universality classes are not included as particular 
cases but only as the limiting cases for i? — )• and i? — )• -|-oo, respectively. 
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For Hamiltonians of Type II, namely n 



n_|_, we have: 



(d) The masses are given by equation (4.53) for k odd and by 



mN,k 



(4.60) 



for k even, where A( J) is as in (4.52). This case clearly belongs to the RP^ ^ 
universality class. 

4.1.5 Interpretation 

We want now to interpret these results in another framework. In one dimen- 
sion a continuum field theory is simply a continuous-time Markov process on the 
target manifold. Now, the generator of a continuous-time Markov process is the 
convex combination of a diffusion part (a second-order elliptic differential operator) 
and a jump part (a positive kernel)."^® Physically, this means that the "particle" 
diffuses for a while according to the specified differential operator, and then, at ex- 
ponentially distributed random times, jumps according to the specified probability 
kernel. On the sphere S^~^ for iV > 3, the only S'0(iV)-invariant second-order 
elliptic operator is the Laplace-Beltrami operator (and multiples thereof); thus, the 
only S'0(iV)-invariant diffusion on S^~^ is standard Brownian motion (with an ar- 
bitrary coefficient, corresponding to a rescaling of time). "^^ On the other hand, there 
is an infinite- dimensional family of possible S'0(iV)-invariant jump kernels K: in- 
deed, one can specify an arbitrary probability distribution of jump angles 6 G [0,7r] 
(and S'0(iV)-invariance then determines K uniquely, for iV > 3). Each one of 
these quantum Hamiltonians H = aC -\- K {a > 0) defines a legitimate continuum 
cr-model. 

Moreover, for each such quantum Hamiltonian H and each t > 0, the integral 
kernel e~^^{(T, cr') is a smooth 0(iV)-invariant function of cr and cr' (and thus a 
function of cr ■ cr'); it can therefore be realized as e^'^""'"" ^ for a suitable smooth 
potential Vt- Thus, by taking some sequence t J, 0, we see that each continuum 
cr-model can be realized as a continuum limit of lattice cr-models (i.e. discrete-time 
0(iV)-invariant random walks on S^~^ ) each of which has a smooth step distribution 



We can now interpret formula (4.46): the continuum limit of this theory is a 
Markov process on S^~^ which contains a jump part with jump angles 9i = arccos ti. 
The coefficients Ai are related to the probability distribution of the jump angles. 
The typical configuration here, for large J, consists of ordered domains where cr^ ■ 
cr^^i 1 separated by links where a jump occurs, that is where cr^ ■ cr^^i ti. 
Notice that these jumps must be sufficiently rare, otherwise they destroy the order 
and thus no criticality appears (this occurs in case 4 of our classification, i.e. when 
ai = a), but not too rare, otherwise they are unable to change the critical behavior 
of the system (this is the case when ai > fi). Jumps of vr (which are simply spin flips) 

^^See [34], Example 1.2.1 (p. 6), Theorem 2.2.1 (p. 48) and Theorem 2.2.2 (p. 51). 

^^This is true also for # = 2 if one demands 0(#)-invariance and not just 5'0(#)-invariance. 



gV(<T.<T')_ 
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play a special role: the spin-A; correlations for k even are insensitive to spin flips, and 
thus they remain critical irrespective of the frequency of such flips. In particular, 
for theories of type II these spin flips are infinitely rapid, and the continuum limit 
is best considered Markov process on 

4.2 Finite-Size-Scaling Limit 

4.2.1 Generalities on the Finite-Size-Scaling Limit 

We want now to discuss the flnite-size-scaling limit for theories of types I and 
II (see Section 4.1.3 for deflnitions). This limit is given by L — )• oo, J — )• +oo 
[hence ( J; oo) — )• oo, where i^^^ denotes any one of the correlation lengths i^^'^''^ 

or i^^Y^ introduced earlier] in such a way that i^^^ (J; oo) /L remains flxed.'*° We 
therefore deflne the scaling variables 

fT r \ — ^N,k^ (J'l /a a-\\ 

Zk = Zk(J;L) = — — . (4.61) 

When considering correlation functions in x-space, we also scale x, i.e. we will 
consider x = IcL with < x < 1 flxed. The corresponding correlations represent 
the correlations of a continuum theory in a periodic box of width 1. 

Everywhere in this section, for theories of Type II or those theories of Type I 
belonging to the RP^~^ universality class [see case I(b2) in Section 4.1.3, and case 
(b) in Section 4.1.4 "^^J, k must of course be even, and in all the formulae below only 
even values of / and m are to be included in the sums. 

As we have seen in Section 4.1.1, the theory displays critical behavior only if the 
masses mjv,fc go to zero in the limit J — )• +oo. Therefore, in this section, we assume 
that 

VN,k{J) = l-aN,kHJ) + o{A{J)) , (4.62) 

where A(J) is a non-universal scale factor (assumed strictly positive) which goes 
to zero as J — )• +00. The quantities ajv,fc characterize the universality class of the 
theory and are completely deflned modulo an overall constant which can be absorbed 
into A( J). Since ujv,fc (J) < 1, we have ajv,fc > 0. For the iV-vector universality class, 
the coefficient ajv,fc can be simply deflned by 

ajv,fc = ^N,k , (4.63) 

■^"^ Similarly to what we did in Section 4.1.1, we could consider a sequence (•)'-"-' of finite- volume 
lattice models with linear lattice sizes ^ cxD. A finit e-size-scaling limit (= finite-volume con- 
tinuum limit) yielding a continuum box of side (0 < < cxd) is defined by rescaling 
lengths by factors S(") = _l(")/_l(<=°"*) cxd) and rescaling field strengths by factors Cjy^ such 
that the spin-A; two-point functions have well-defined limits. Without loss of generality we can set 

J^(cont) _ ^ 

'^^In these two cases, the masses of the even and odd sectors go to zero with different rates. As 
explained in the remark at the end of Section 4.1.1, A(J) is chosen to be of the order of the smallest 
mass of the theory (that of the even sector) and it can be seen that for k odd t;jv,fc ( J) goes to zero 
exponentially in the finite-size-scaling limit. Therefore, the odd sector does not contribute to the 
finite-size-scaling functions, just as for case (a) of the RP^~^ universality class. 
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where Ajv,fc are the eigenvalues of the Laplace-Behrami operator on the sphere. For 
the RP^~^ universahty class the same holds for even k. For the other universality 
classes, which include jump processes, the coefficients ajv,fc can be easily derived 
from (4.43) for theories of Type I and its analogue (4.51) for theories of Type II. 
Now, from (4.61)-(4.62) and (3.24)/(3.26), it follows that, for large J, we have 

zk ^ ll[m,kLK{J)] . (4.64) 

Therefore, instead of considering the limit L, J — )• oo Zk fixed, we will equiva- 
lently consider the more convenient limit at L K{ J) = 7 fixed, since the parameter 
7 will appear naturally in our formulae. Let us then define the variables 

= -,( J; L) = ^ = . (4.65) 

(To leading order"*^ we have z^ ~ The distinction between z^ and will become 
relevant only when we consider corrections to finite-size scaling.) Our approach 
will be to compute various quantities as a function of 7, and then use (4.65) to 
re-express everything as a function of Ji or "2- The reason for this last step is 
that functions of 7 are universal only modulo a scale factor [corresponding to the 
arbitrariness of A( J)], while functions of physical continuum quantities (such as the 
Ik) are universal tout court. 

4.2.2 Computation of the Finite-Size-Scaling Functions 

We want to compute the following finite-size-scaling functions:'*"^ 

ZS,°)(7) = lim ZNiJ;L)/FN,oiJf (4-66) 



L, J ^ 00 
7 fixed 



Gj^),(x,7) = lim GnA^L,J;L) (4.67) 

L, J CO 



7 fixed 



L, J ^ 00 
7 fixed 



XN'kiV = . Ijm (4.68) 



42.^)(o)(^) ^ lin, ^^-^ Y (4-69) 

Lj 



J ^ CO 
7 fixed 
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More precisely we have 



lim _ , ^' = Jim CnY '^N,k MJ) = 1 



■*^The superscript ("^^ indicates "leading order" . The first corrections to these finite-size-scaling 
functions will be computed in Section 4.3. 
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Concerning XN,k and it is often convenient to look at the ratios 

R.;nAJ:L) ^ MAIiR (4.70) 



Xjv,fc(J; oo) 



Rc;nAJ;L) ^ . (4-71) 



Note that by (3.22)/(3.23) we have 



Xjv,fc(J;^) ^ T—- = 2Lzk (4.72) 

hence the ratios have weU-behaved finite-size-scahng hmits: 

R%A-f) = hm R^,^,,{J-L) = \~a^,,-fx%]k{l) (4-74) 

7 fixed 



i?[°U7) = lim R,,nAJ;L) = ~a^,kl&^'\^) (4.75) 

L, J ^ oo 
7 fixed 



The computation of the finite-size-scahng functions (4.66)-(4.69) is straightfor- 
ward. In the hmit L, J — )• oo with 7 fixed, we have 

VN,k{Jt = VN,k{jy'^^'''' ~ exp I ^^^jy log [1 - Sjv^fc A(J)]| exp(-7ajv,fc) . 

(4.76) 

Inserting this limit in the exact expressions (3.13)-(3.15) from Section 3.1, we obtain 

00 ^ 

ZPil) = J2^N,ie-™ (4.77) 

1 00 r'2 ^ 

G^Z(-rf) = -(5)— E e--.' e--^-^ (4.78) 

Zyi-f) l,m=0 ^N,k 

{00 r 1 ^ 1/2 
E [BN;k,l,mil) / K;l,m\ 
'-^^^^ } (4.80) 
E BN;k,l,mi7) 
Lm=0 
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AjV;/,m = aN^m — aN,i (4-81) 
Bn;,,Ui) ^ Ki,m . , ^^""^ % (4.82) 

Let us notice that (4.77) can be rewritten as 

4°) (7) = Tr exp(-7:ff), (4.83) 

where H is the operator that generates the continuous-time Markov process corre- 
sponding to that universahty class. For Type-I theories (except the case belonging 
to the RP^~^ universality class), the trace is taken in the space L'^{S^~^); while 
for Type-II theories (and for the case of Type I which falls in the RP^~^ univer- 
sality class), the trace is taken in the space L^{RP^~^), which is isomorphic to 
L'^{S^~^)even = 0fcLo, k even EN,k and cousists of the tven functions on S^~^ . Phys- 
ically, (4.83) expresses the fact that the finite-size-scaling limit corresponds to the 
continuum theory in a finite periodic box. 

Notice that since the coefficients ajv,fc are uniquely defined by the universality 
class of the theory, modulo a A;-independent rescaling [which depends on the explicit 
definition of the scaling factor A(J) but does not affect the products ajv,fc7], these 
finite-size-scaling functions are universal modulo a rescaling of 7. 

4.2.3 An Interesting Family of Universality Classes 

Let us examine in more detail the finite-size-scaling curves for R^-]\j ^k{J ] L) . In 
particular, we want to study their dependence on the different universality classes 
described in Section 4.1. As can be seen from the explicit expression of -R^°jv kW)i 
finite-size-scaling curve is determined completely by {Sjv,/}. Therefore, we consider 
a family of universality classes parametrized by a continuous variable 5, with Sjv,/ 
given by 

r Ajv,; for / even 

"^'' = |A^., + i? for /odd ^^'^^^ 

(This family of universality classes was found in part (c) of the example in Section 
4.1.4, and will be also found for the two-parameter Hamiltonians treated in Section 
5.) We can get the iV-vector universality class by choosing 5 = 0, and the RP^~^ 
universality class by taking the limit B ^ 00. 



■^■^As mentioned in footnote 21 above, the expression (3.15) for XN,k requires some exegesis when- 
ever iijv,; = VN,m', and correspondingly (4.79)/(4.80) require exegesis whenever ajv,; = o.N,m- In such 
cases the combination e~'''""''/Ajv; ;,m, which occurs in (4.79) and in the numerator of (4.80), should 
be interpreted as ■je~'^'^'^'' /2. This can be seen by going back to (3.14); it can also be obtained by the 
"quick-and-dirty" method of symmetrizing in / and m, treating the ajv,; as if they were independent 
variables, and using I'Hopital's rule. Note that in the vector and RP^~^ universality classes this 
problem occurs only when / = m (and hence k is even). However, in the more general case (4.84), 
for certain values of B one may have ajv,; = o.N,m for / 7^ m (but only where / — m and k are odd). 
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Let us first look at the isovector sector (k = 1). In Figure 1 we plot Rl^.jsj -^^{ j) 
for various values of the parameter 5, using as an example iV = 4. The graphs 
are drawn not as functions of 7, but as functions of the more "natural" variables 
Jk = l/(ajv,fc7) defined in equation (4.65). In Figure 1(a) we plot versus Ji, while 
in Figure 1(b) we plot versus "2! different aspects of the behavior can be observed 
in these two plots. 

A few interesting features that can be seen in the graphs for iV = 4, and that 
can in fact be proven easily for arbitrary iV, are: 

(i) In the limit 7 — )• (i.e. Ji, "2 00), we have lim i^^^jy ^(t) = (for finite B). 
More precisely, an expansion for small 7 of -R^°jvfc(7i ^) arbitrary k gives 

- ^ + 0(1) . (4.85, 

independent of B. This behavior is observed in Figure 1(a), where the dashed 
curve represents (4.85). 

(ii) For < -B < 2 the curve is decreasing at small Ji (or J2), while for B > 2 it 
is increasing: this can be seen from a large-7 expansion of i?^°]Yi(7)- 

(iii) lim i?[°]y-^(7; B) = 1 for all fixed 7 > (i.e. all fixed J2 < 00). This behavior 
is observed in Figure 1(b). 



Let us next look at the isotensor sector (k = 2). In Figure 2 we plot the ratio 
^x-^N 2(1) ^ function of "2 for three different values of the parameter B, for the 
case iV = 4. Figures 2(a) and (b) show the same curves, but emphasizing different 
ranges of the variable "2- A few features that can be seen in the graphs for iV = 4, 
and that can also be checked from the explicit formulae for general iV, are the 
following: 

(i) The curves are monotonically decreasing functions of the family parameter B 
for each fixed value of the abscissa "2- [We can write 

R%Ar, B) = R%,,{r, 0) - [R%,,{r, o) - R%,,{r, ^)] 

where Z^^ (or 4°^^) is Z^"^ with the sum restricted to even (or odd) /, and 
all the Zjv's are evaluated at i? = 0. This proves the monotonicity in B for 
fixed 7.] 
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FSS function for ^N.k- N=4, 1^=1 

— I r — I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 




0.0 



1 2 

Zl 



FSS function for Xn.ic^ N=4, k=l 



K 0.5 




Figure 1: Graph of the ratio -R^°jvi ^ function of (a) Ji and (b) "2, for the 
case iV = 4, for the family (4.84) of universahty classes. In (a) the lowest curve 
corresponds to i? = 0, which is the iV-vector universality class; the highest curve 
is i? = 20; the third solid curve is the limit B — )• +00; and the dashed curve is 
the asymptotic behavior (4.85). In (b), the lowest curve is i? = 0; the next three 
curves are B = 2^ B = 8 and B = 20, respectively; and the straight line is the limit 
B +00. 
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FSS function for Xn.r^ N = 4, k = Z 




Q Q I \ \ \ \ I \ \ \ \ I \ \ \ \ I \ \ \ \ I 

0.0 0.2 0.4 0.6 0.8 



FSS function for XN.k^ N=4, k = 2 




Q Q I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 

12 3 4 



Figure 2: Graph of the ratio -R^.jv2 ^ function of "2, for the case iV = 4, for 
the family (4.84) of universahty classes. In (a) the highest curve is i? = 0, and 
corresponds to the iV-vector universality class; the lowest curve is i? = 00, and 
corresponds to the the RP^~^ universality class; the curve in-between is i? = 1. In 
(b) the higher curve is obtained for B = and the other for B = +00; the dashed 
curve is the asymptotic behavior (4.85). 
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(ii) The curves coincide exponentially rapidly for large Z2 (i.e. small 7). Indeed, 
for k even, the behavior for small 7 is (see Appendix B.6) 

41. (7; B) = ^ PMl) [1 + 0(e--V^-)] , (4.87) 

where Pjv,fc(7) is a polynomial independent of B. (More precisely, in Appendix 
B we shall prove this behavior only for A; = 2, but we conjecture that it 
holds for all even k.) This is why the dependence on B disappears in Figure 
2(b) long before the curves show the asymptotic behavior (4.85). Physically, 
the behavior (4.87) reflects the fact that the universality classes (4.84) are 
equivalent at all orders of perturbation theory; the i?-dependence is a wholly 
nonperturbative effect. A similar situation occurs in the iwo -dimensional a- 
models [5,6,7]. 

(iii) The curve for the RP^~^ case is not monotonically decreasing as a function 
of J2, but is slightly increasing for small J2. (In fact, an expansion for large 7 
shows that the function is increasing at small J2 for all values of B.) 

Finally, we show the kind of flnite-size-scaling plot that one usually considers 
in Monte Carlo simulations: here the inflnite- volume correlation lengths are not 
known, so instead of Zk or we would use the variable 

X, = x,{J-L) ^ Wl^lilZ. (4.88) 

Moreover, in this case we cannot compare lattice size L with 00 [as requested in 
(4.70)]; rather, we must compare L with (for example) 2L [11,13]. In Figure 3(a,b) 
we show the analogues of Figure l(a,b): that is, we plot the flnite-size-scaling 
curves for the ratio Xn,i{J] L)/xn,i{J] 2-L) as a function of Xi and X2^ respectively, 
for various values of B. The FSS curve for this ratio is given by XN\il) / [^XN\i^l)] 
plotted parametrically versus 

Xk ^ lim x,{J;L) = d'f(°^(7) • (4.89) 

L, J ^ 00 
7 fixed 

In Figure 4 we show the analogous plot for the isotensor sector. 

It is interesting to compare the curves in Figures 3(a) and 4 with those coming 
from a Monte Carlo study of a similar family of universality classes in two dimen- 
sions [5,6,7]. The curves are qualitatively very similar, although of course they are 
quantitatively different. 

Now let us compare the flnite-size-scaling curves to the explicit solution at flnite 
L and J. We show in Figure 5 the graph of the flnite-size-scaling function of the 
spin-1 susceptibility for the iV = 4 iV- vector universality class [namely the lowest 
curve in Figure 1(a)] together with some points calculated from the exact expression 
of Ry,-N^i(J; L), for several values of L. More precisely, we plot: 
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0.0 ^ 
0.0 



0.5 



1.0 



1.5 



L/2L FSS function for Xn.r^ N = 4, k=l 




Q Q I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 

0.0 0.2 0.4 0.6 



Figure 3: Graph of the ratio XN\il) / [^XN\i^l)] ^ function of (a) Xi and (b) X2, 
for the case iV = 4, for the family (4.84) of universahty classes. 
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L/2L FSS function for Xn.i,: ^=4, k = 2 
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XNk(zj^): N = 4 N-vector model, k=l 
I I I I I I I I I I I I I I I I I I I I I 




Q Q I I \ I \ \ \ I I \ \ \ \ I I \ \ I \ \ I 

0.0 0.5 1.0 1.5 2.0 

Figure 5: Graph of Ry,-N^i(J; L) as a function of Ji = (Sjv,i7)~^ for the one- 
dimensional iV = 4 iV-vector modeL Symbols indicate: L = 4 (*), 8 (+), 16 
(x), 32 (□). The corresponding finite-size-scaling function is also plotted. 



• Points: exact values of Ry,-N^i(J; L) plotted versus Ji = l/(ajv,ii^A( J)) for 
the iV-vector model (3.2), using the formulae from Section 3.1 and ujv,fc(J) 
given by the Bessel functions in (2.32). Here A(J) = 1/(2 J) and ajv,fc = Ajv,fc- 

• Curve: the finite-size-scaling function i?^°]Yi(7) for the iV-vector universality 
class, function of 2;i = (ajv,i7) ^. 

Notice that for small values of L there are significant corrections to finite-size scaling. 
These corrections will be discussed in the next section, and we will show that they 
are of order 1/X, or equivalently, of order A(J). 

Let us now return to the original (and most natural) scaling variables Zk = 
^nY\ J'i ^)/ L. The finite-size-scaling curves are of course the same, since Zk and Ik 
coincide at leading order. However, the meaning of the points in the plot is different, 
and the corrections to finite-size scaling may differ in the two variables. In Figure 6 
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XNk(zj^): N = 4 N-vector model, k=l 
I I I I I I I I I I I I I I I I I I I I I 




Q Q I I \ I \ \ \ I I \ \ \ \ I I \ \ I \ \ I 

0.0 0.5 1.0 1.5 2.0 

Figure 6: Graph of R^-n^i( J; L) as a function of Zi = ^^^^^\ j]oo)lL for the one- 
dimensional N = A iV-vector modeL Symbols indicate: L = A (*), 8 (+), 16 (x), 
32 (□). The corresponding finite-size-scaling function is also plotted. 



we show the same data points as in Figure 5, but plotted versus Zi instead of Ji. In 
Figure 7 we make the analogous plot for the "L/2L" FSS function plotted versus 
Xi (which is a close relative of Zi). Clearly, the corrections to finite-size scaling are 
considerably smaller if we use variables Zk or Xk rather than the variables Ik- In 
Section 4.3 we will show that the 1/L corrections vanish in the variables Zk or Xk] 
the leading correction appears to be of order 

The fact that the plot in terms of Zk or Xk shows better agreement with the 
finite-size-scaling curve than the plot in terms of can be interpreted as a mani- 
festation of the difference between "scaling" and "asymptotic scaling". As used by 
lattice quantum field theorists, these terms mean the following (see e.g. [7]): "Scal- 
ing" denotes the convergence to the continuum limit for dimensionless ratios of 
long-distance observables and for the relations between such observables. "Asymp- 
totic scaling", by contrast, denotes the convergence to the asymptotic predictions 
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XN_k(xi^): N = 4 N-vector model, k=l 

I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I 




^< — i — >^-fi-^^ — \ — 



Q Q I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ I 

0.0 0.5 1.0 1.5 2.0 

Figure 7: Graph of Xn,i{ J', L)/ Xn,i{ J', ^L) as a function of Xi = Cn'i^\ J', L)/L for 
the one-dimensional iV-vector model with iV = 4. Symbols indicate: L = 4 (*), 8 
(+), 16 (x), 32 (□). The corresponding finite-size-scaling function is also plotted. 
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(e.g. as J — )• oo) for the relation between long-distance observables (such as x 
or 1^ or combinations thereof) and the "bare" parameters in the Hamiltonian (i.e. 
J)."*^ Clearly, asymptotic scaling (to a given degree of accuracy) implies scaling (to 
the same degree of accuracy), but not conversely; otherwise put, the corrections 
to asymptotic scaling may be much larger than the corrections to scaling. Now, 
Rx\N,k{J]L) = XN,kiJ'^ L)/ XN,kiJ'^ and Zk{J;L) = ^^J^(J;oo)/L are examples 
of dimensionless ratios of long-distance observables, and the rapid convergence to 
the continuum limit in Figure 6 is an example of rapid "scaling". On the other 
hand, 'Zk{J]L) = l/[ajv,fc A( J)] is an example of a bare parameter (by virtue of 
its explicit dependence on J), and the slower convergence to the continuum limit 
in Figure 5 is an example of the less-rapid onset of "asymptotic scaling". In this 
model the corrections to asymptotic scaling are of order while the corrections 
to scaling appear to be of order 

4.3 Corrections to Finite-Size Scaling 

In this section we want to compute the corrections to the finite-size-scaling 
functions. We assume a large- J expansion of the form"*^ 

VN,k{J) = 1- aN,k HJ) + bN,kKorriJ) + o{AcorriJ)) , (4.90) 

where Acorr{J)/ A{J) goes to zero for J — )• -|-oo. 

In the limit J — )• -|-oo, L — )• oo at LK{ J) = 7 fixed we have 



1 I U ^corr\J) 7-2 K / j\ 

1 + 7 0iv,fc — 777^ 7«JVfcA(-') 




VN,k{J)^ = exp(-7ajv,fc) 



A(J) 

= exp(-7SAr,fc) [1 -76jv,fcA(J) + o(A(J))] , (4.91) 

where A(J) is the more slowly decreasing of Acorr{J)/ A{J) and A(J), and — 6jv,fc 
is the corresponding coefficient. (Of course, if A(J) and Kcorr{J) I A.{J ) are of the 
same order, then — 6jv,fc is given by the sum of the two coefficients.) Plugging this 




■*®In place of the bare parameters, one may alternatively use short- dt stance quantities such as the 
energy E , inasmuch as they play a similar physical role. 

■*®In Section 4.1.2 we computed aN,k and &jv,fc for two simple cases of Hamiltonians h{t): 

(i) t = 1 is the only absolute maximum and > 0. (This generalizes the #-vector model.) 

(ii) t = ±1 are the only absolute maxima, h{t) is an even function and /i'(l) > 0. (This is the 
symmetric subcase of what we called in Section 4.1.2 the "second simple case": see footnote 
33 and the text following it. It generalizes the RP^~^ model.) As discussed before, we look 
only at the even-A; sector in this case. 

In both simple cases we obtain the same coefficients aAr^, &jv,fc [cf- (4.19)-(4.21)]; and we have 

A(J) ~ 1/J, Aeorr(J) ~ 1 1 P ■ 
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expression in (3.13)-(3.15) we immediately obtain the corrections to the finite-size- 
scaling functions. Notice that the procedure is straightforward, and by adding more 
terms in the expansion (4.90) we can compute the corrections to any arbitrary order. 
For example for the susceptibility we obtain 



XN,AJ;L) = L Ix^Ut) + xIvUt) A(^) + oiA{J) 



(1) 



(4.92) 



where 



XnI (t) 



E 



c 



N; k,Lm ^ 



7 Zjj^ (7) l,m=0 ^N,k AjV; l^rn 



A 



7 b 



NA 



N: Lm 



+ 



7 



(7) 



Zat' \ T ) n=0 

In the same way, the ratio i?^;jv,fc(<^; L) defined in (4.70) is given by 

Rx.,nAJ:L) = R%,kii) + MJ)R%,,ii) + o(a(j) 

with 



^■N,k (7) 



7 ~ 



XN^k (7) + ^XN^k (7) 



where we have used (4.74) and the expansion 



Xjv,fc(J;oo) 



aN,k M/) 



JN,k 



(lN,k 



A(J) + o A(J) 



(4.93) 

(4.94) 
(4.95) 

(4.96) 



These formulae simplify considerably in case 6jv,fc has the simple structure 



bN,k — (lN,k Cn 

for some coefficient cjv- In this case we obtain the simple formulae 

XN^k (7) 



d 



R 



(1) 

x;N,k 



(7) 



cjv7^X]v,fc(7) 

7 ~ d ^ (0) / ^ 

2 «iv,fc CJV ^ [7 XN,k (7) 



(4.97) 

(4.98) 
(4.99) 



In particular, the two simple cases mentioned above (see footnote 46) satisfy this 
requirement with 

CJV = (iV + l)r + 1 . (4.100) 

Indeed, we have here Acorr{J) = HJf, so that A(J) = A(J) = 1 1 [2h' {!).]] and 
^JV,fc = k ~ ^N,k', the claim then follows from (4.19)-(4.21). Note that -R^^.jvi (7) 
here depends on r only through the global factor (N -|- l)r -|- 1. 
We now study in more detail the classes satisfying (4.97). 

In Figure 8 we show the correction to finite-size scaling for the spin-1 suscepti- 
bility for the N = 4,8 iV-vector universality class (cjv given by (4.100) with r = 0). 
We plot function of Zi [defined in (4.65)] 
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Figure 8: Corrections [i?^;jv,i { J] L) — -R^.jvi (t)] to the finite-size-scaling function 
of Ry,-N^i(J; L) for the one-dimensional iV-vector model as a function of Ji. The 
higher and lower curves correspond to iV = 4 and iV = 8, respectively. Symbols 
indicate L = 4 (*), 8 (+), 16 (x), 32 (□). The function (4.99) (for these two values 
of N) is also plotted. 



• Points: the difference 

[R^.,nAJ:L) -R%,,{i)]L (4.101) 
for the iV-vector model (with iV = 4, 8 and for different lattice sizes) and 

• Curves: the corresponding limiting curves -R^^.jvi (t) given by (4.99)/(4.100) 
with r = 0. 

Now we can compare the corrections shown in Figure 8 to the corrections for the 

(-U.\ 

case in which we choose Zi = i^j^i (Jj °°) / L as the variable in the abscissa. As was 
done in the previous section, we must replace 7 in formula (4.94) by its expansion 



44 



in terms of Zk to the desired order, given by 



7 



We get 



R 



(lN,k Zk 



{J-L) 



bN,k — 



A(J) + o A(J) 



R 



(0) 

x;N,i 



+ o A(J) 



(4.102) 



(4.103) 



That is, there are no corrections at order A(J) in this case; in other words, the 
corrections of order 1/L found for Figure 8 are not present here! EmpiricaUy it 
appears that the leading corrections are in fact of order see Figure 9. (The 

hmiting curve shown was evaluated numerically by taking a very large value of L.) 
The same holds for the corrections to finite-size scaling written in terms of the 
variable Xk defined in the previous section, as shown in Figure 10. As discussed 
in Section 4.2.3, this is a manifestation of the difference between "scaling" and 
"asymptotic scaling". 



5 A Two-Parameter Family of Hamiltonians 

In the previous section we have investigated the continuum limits arising from 
a one-parameter family of interactions. One might imagine that, by considering 
many-parameter families of Hamiltonians and taking appropriate trajectories in 
the multi-parameter space, one could find additional continuum limits. We have 
investigated this problem for a two-parameter family of interactions given by 

h{(T ■ t) = Jv hv{(T ■ t) + Jt hxicr ■ t) . (5-1) 

We will not study the problem for generic hy and /ix, but will restrict our discussion 
to the case in which hy is an odd function and has a unique maximum at 1 while 
hx is an even function and has maxima at ±1. Moreover, we will assume hy(l) > 
and h'j'(l) > and we will consider only the case Jy, Jt > 0. This generalizes the 
mixed isovector/isotensor model 

h{a-T) = Jv(T-T + ^{(T-rf (5.2) 

studied in [5,6,7]. 

We want now to find the critical points of these theories. Since in dimension 
d = 1 no phase transition can occur for finite values of the couplings, we must 
investigate the limit in which at least one of the two couplings tends to infinity. It 
is trivial to see that in the limit Jj — )• -|-oo with Jy fixed and finite, one recovers the 
JlpN-i universality class; while in the limit Jy +oo with Jj fixed and finite, one 
reobtains the iV-vector universality class. It therefore remains only to investigate 
the case in which both Jy and Jt go to infinity. 
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Figure 9: Corrections [i?^;jv,i {J] L) — R-^.j^ i (7)] to the finite-size-scaling function 
of Ry,-N^i(J; L) for the one-dimensional iV-vector model as a function of Zi. The 
higher and lower curves correspond to iV = 4 and iV = 8, respectively. Symbols 
indicate: L = 4 (*), 8 (+), 16 (x), 32 (□). The corresponding limiting curve 
(numerically evaluated for a large value of L) is shown for both cases. 
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Corrections to FSS of Xni(^i) 

0.3 I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r 




0.0 0.5 1.0 1.5 2.0 



Figure 10: Graph of { [xjv,i (J; L) /xn,i {J; 2i:)]-[x]v,i(7)/Xjv,i(2 7)]}^^ as a function 
of Xi for the one-dimensional iV = 4 iV-vector modeL Symbols indicate: L = 4 (*), 
8 (+), 16 (x), 32 (□). The corresponding limiting curve (numerically evaluated for 
a large value of L) is also shown. 
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(a) Let us first consider trajectories sucli tliat Jfj Jy — > as Jy, Jx — > +00. In 
tliis case, from (2.27), we get 



+ O 



1 J7 



''N,k 



(5.3) 



wliere 



fni-Jv-, Jt) 



Mi) 



:27r/i'(l)l^/' V2 
1 



l-JV/2 



1 iV^ - 1 
\N,k + - (iV - 1) (iV - 3) ^ ry 



(5.4) 



(5.5) 



2hy{l) 

wliere Ajv,fc are tlie eigenvalues of the Laplace-Beltrami operator on the sphere and 

rv = hy{l)/h'y{l) . (5.6) 

Thus, in this limit hx is an irrelevant perturbation, and we get the iV-vector uni- 
versality class. 

(b) Next let us consider trajectories such that Jt/Jv = with < a < 00. In 
this case we can rewrite (5.1) as 

(5.7) 



h(a ■ t) = Jv \ hv (cr ■ t) + ahx (<t • r) . 



This is a one-parameter family of interactions with Hamiltonian h which has a 
unique maximum at t = 1. Thus also in this case we get the iV-vector universality 
class. 

(c) Finally, let us consider trajectories such that Jx/Jy — > 00. We get from 
(2.27) 

FnAJv.Jt) = /iv(Jy,JT) I 1-^ + 0(^-1,^^ + 



+ (-l)^exp -2Jvhv{l) 



1 + 



(5.8) 



where /jv(Jy, Jt) is defined in (5.4) 

1 



and 

It follows that 



2/iy(l) 



1 iV^ - 1 
Ajv,fc + - (iV - 1) (iV - 3) ^ rr 



r-T = /i;(i)//i;(i) . 



VN,k{Jv. Jt) ~ 1 - ^ - [l-(-l)'l exp[-2Jy/.y(l)] 



(5.9) 
(5.10) 
(5.11) 
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where 



A 



ojvfc = ~, ' — . (5.12) 

To go further we must distinguish three different cases according to the relative size 
of the two correction terms in (5.11), i.e. according to the behavior of the product 
JTexp[— 2Jy/iy(l)]. 

(i) Let us first consider trajectories for which Jx exp[— 2 Jy/iy(l)] goes to zero. In 
this case the exponential term in (5.11) goes to zero faster than the 1/ Jt term 
and can thus be dropped. We reobtain in this way the iV-vector universality 
class. 

(ii) In the opposite case, when Jx exp[— 2 Jy/iy(l)] — )• +oo, the leading behavior 
is given by 



VN,k ( Jv, Jt) 



so that 



mN,k 



{1 — (iN,k I Jt for k even 

1-2 exp[-2Jy/iy(l)] for k odd 

o^N,k I Jt for k even 

2 exp[— 2 Jy/iy(l)] for k odd 



Thus for all k odd 

T^N^k 



mN,2 



(lN,2 



exp[-2Jvhv(l)] 



mN,k 2J 

so that these limits belong to the RP^~^ universality class, 
(iii) Finally let us suppose 



(5.13) 



(5.14) 



(5.15) 



Jt exp[—2Jvhv{l)] 

where i? is a constant. In this case we get 

C 1 — Sjv^fc/ Jt 



B 



VN,k ( Jv, Jt) 



1 



(lN,k + 



au' 



d 



mN,k 



4hj.{l) 



Ahjil) 

for k even 
^ Jt for k odd 

for k even 



(5.16) 



2/iy(l) Jt I ^N,k + B for k odd 



(5.17) 



(5.18) 



Thus, for < -B < oo we find again the intermediate universality class (4.58) 

— interpolating between the iV-vector and the RP^~^ universality classes 

— which appeared already for one-parameter Hamiltonians with maxima at 
t = ±1. 
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Let us notice that in this last case the even-spin correlation functions are equal 
to those of the iV-vector model, while the odd-spin ones are a product of an Ising 
correlation and the corresponding iV-vector correlation. This family of theories 
is parametrized by i? (0 < -B < oo) and all the limiting mass ratios TZ^^k are 
determined in terms of B as in equation (4.50). Equivalently, we can choose any 
one of these ratios (with k odd) to characterize the universality class; for instance 
we can use the ratio 

In the continuum limit we have 

Thus each theory is labeled by the ratio 7^jv,i , which can assume any value from to 
2N/(N — 1). Notice one special feature of d = 1: the maximum value of mjv,2/"^JV,i 
is not 2 but rather is larger. This is due to the fact that in (spacetime) dimension 
d = 1 scattering states cannot exist, so the usual inequality mjv,2 < ^"m-N,! does not 
apply. 



A Properties of Hyperspherical Harmonics 
A.l Calculation of AfN,k = dim Ej^^k 

Let us begin by computing the dimension of the linear space -Ejv,fc consisting of 
the completely symmetric and traceless tensors of rank k over R^. This can be 
done by computing the dimension of the space of all completely symmetric tensors 
of rank A;, and then subtracting from it the number of independent trace conditions 
that have to be imposed to ensure the tracelessness of these tensors. The number of 
linearly independent symmetric tensors is given by (^"^^"^^ (the number of ways of 

placing k prisoners in N cells), and the number of traces is given by (the 
same binomial as before but considering only k — 2 indices; of course this simply 
vanishes if A; < 2). Therefore we obtain 

1- ^ /iV + A;-l\ /iV + A;-3\ , 

^N,k = dim EN,k = , - , o (A-la) 



(A. lb) 



T{N + k) T{N + k-2) 

k\r{N) ~ {k-2)\r{N) 

N + 2k-2 T{N + k-2) 
k\ T{N - 1) 



(A.lc) 



[with the interpretation (—2)! = ( — 1)! = oo in (A. lb)]. This proves formula (2.2). 

We shall take (A.lb)/(A.lc) as the definition of A/jv,fc for N noninteger. (By 
contrast, we shall always consider k to be an integer > 0.) Note that for each fixed 



50 



integer A; > 0, A/jv,fc is a polynomial of degree k in N; in particular, it is well-defined 
and finite for all real N. Note also that for each fixed N (not necessarily integer), 
we have A/V.fc ~ 2A;^-Vr(iV - 1) as A; 
The simple identity 



oo. 



AT, 



JV+2,fc-l 



A;(iV + A;-2) 
N{N - 1) 



AT, 



N,k 



(A.2) 



will play an important role in Appendix B. 

Finally, for integer N > 3 we have the following formula: 



AT, 



' (JV-4)/2 

n 

r=0 



k + 



N-2 



N,k 



{N -2)\ 



X < 



, (N-5)/2 r 

k + ^y n 



JV-2 



r + i 



for N even > 4 



for iV odd > 3 



(A.3) 

In particular, for even (resp. odd) iV > 3, A/jv,fc is an even (resp. odd) polynomial in 
the "shifted index" k + ^. Note, however, that for iV = 2, ^^N, k is not a polynomial 
in A;: for A; > 1 we have A/jv,fc = 2, in agreement with (A.3), but A/jv,o = 1 7^ 2. 



A.2 Some Basic Formulae 

Let us now compute the integral of a product of an even number of <t's (an 
odd number gives trivially zero). Let us introduce, for an arbitrary vector A„, the 
quantity 

h{A) = J dflia) (A-af . (A.4) 

As dQ (<t) is rotationally invariant, we have Ik (RA) = Ik (A) for every R G 0(iV), 
so Ik{A) depends only on Moreover, Ik is manifestly a homogeneous function 
of degree 2 A;. Hence we must have Ik (A) = Jk [A^]^ for some constant Jk- Now, as 
<T^ = 1, we get from (A.4) 

Ik (A) = 2k {2k -I) Ik-i{A) (A.5) 



5A„ dA_ 

A recursion relation for Jk immediately follows 



2k — 1 

= — Jfc_i (A. 6) 

iV + 2A; - 2 ' 



Using Jo = 1 we obtain the general solution 



F A; + l F f 
Jk = 7^ (A.7) 

Fi Ff + A;^ 
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Taking then 2k derivatives with respect to A in (A. 4) we obtain the weU-known 
resuh 

^ ^ 2^r(f + A;) ^ ^ ^ ^ 

where the terms in parentheses correspond to the (2k — 1)!! different pairings of the 
indices. 

Let us now prove the orthogonahty relation (2.10). This is completely equivalent 
to proving that for arbitrary completely symmetric and traceless tensors Tjv,fc and 
Un I we have 



rj^ai^.au JjPl-Pl _ J'j^^^ . If^^k ■ (A-.9) 



To see this, let us first use the definition (2.3) and let us notice that the "Traces" 
terms do not give any contribution due to the tracelessness of Tjv,fc and Un^i. Thus 
the l.h.s. in (A. 9) becomes simply 

I^N,k I^N,i J dflia) a"\.. ...o^' T^;,-"'= 1]^^,^'^ . (A.IO) 

Then let us use (A. 8). The only non- vanishing contributions come from those terms 
which do not contain (5"'"^ or (^'^''^j; such terms exist only if / = A;. In this last case 
there are k\ equivalent contractions and we end up with 



k 



r(f 

^' k\ 



2'=r(f + A; 



We thus obtain the orthogonality relation (2.10) for the K's, provided that they are 
normalized as in (2.4). 

The general formula for the hyperspherical harmonics'*'' can be obtained by using 
the fact that they are completely symmetric and traceless. The complete symme- 
try, together with the needed transformation properties under S'O(iV), implies an 
expansion of the form 

Lfc/2J 

y^^T'i'^) = w E ^iv,^;. ^(i:;r l*^) , (A.12) 

s=0 

where 

Pikisf i^) = ■ ■ ■ S"^'-'"^' . . . a"* + permutations (A. 13) 

and the number of permutations necessary to make P(k;s) completely symmetric is 
(2s — 1)!! (^2*1) • Now we impose the tracelessness. We first note that 

P^;^" (^) = P^:^^ (^) + (iV + 2A; - 23 - 2) P"-- [a) , (A.14) 

■*^This result is obtained in [27]. Note that in their notation P includes all the k\ permutations, 
i.e. it is [{k — 2s)! s\ 2"]"-^ times our P. 
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and therefore we get, from '^aia2^jvfc ("") ~ recursion relation 

AN^k;s-l + {N + 2k-2s-2) AN^k;s = . 

Imposing the normahzation Ajv,fc;o = 1, we find 

{-ly r(f + A;-3-i) 



A 



N,k:. 



Thus we can write 



[k/2\ r(f + k-S-l 

Yn, (-) - f^N, g ^ r(f + ^-1^ "^^^^^^ 



(A.15) 



(A.16) 



(A.17) 



Let us now discuss the relation between the hyperspherical harmonics and the 
Gegenbauer polynomials. From Section 2 we know that Yj^'^^ (<t) is the restriction 
to the unit sphere of a degree-A; harmonic polynomial. Moreover it depends only on 
cr^, so that the polynomial can be written as r^Pk (xi/r) where r = |x|. Requiring 
the polynomial to satisfy Laplace's equation we get for Pk{x) the equation 



X (N-1)-^ +k(N + k-2)Pk 
ax 







The regular solution of this equation'*® is the Gegenbauer polynomial C\ 
The normalization is fixed by the requirement that 



(A.18) 



Y^'l'' (cr) = /ijv,fc(cr"^)'^ + lower-order terms 



We thus get'*^ 

which, using the fact that 



f + ^-1 



N 



v4f-l. 



c 



N/2-1 



(1) 



'N + k-3^ 

, k , 



(A.19) 



(A.20) 



(A.21) 



gives (2.23). 

Note that we could have derived (A.17) by using (2.23) and the expansion of 
the Gegenbauer polynomials [see [19], formula 14 (p. 294)]. 



^^See [31], p. 1031. 

■*^Note that this formula is well-defined in the limit N ^ 2. See footnote 14. 
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A. 3 The Projector onto Symmetric Traceless Tensors 

Using properties 1 to 4 of the projector i'^-"i',fh-fh derive its general 

expression. We start by noting that the most general form satisfying the symmetry 
properties 1 and 2 is 



J-N,k - 1^ ^NMsQ(k;s) ' [A.ZZ) 

where 

Qai.--ak;Pl---Pk — ^0.10.2 ^a2s-i(y.2s ^PlP2 ^/?2s-l/?2s ^a.2s + iP2s + l ^(y-kPk 

+ permutations (A. 23) 

[i.e. there are s S^s among the a's, s among the /3's, and k — 2s connecting the a's 
with the /3's] and the number of permutations necessary to make Q(k;s) completely 
symmetric is given by 



k\ V 1 



(A. 24) 



\s\2'J {k-2s)\ 
Now notice that a consequence of properties 3 and 4 is 

T-ai •••"A: ;/?!•••& vf^l-'-f^k _ ,, T-ai •••"A: ;/?!•••& ^Pi _ -s^ai...ak /A r)t-\ 

^N,k ^N,k — I^N,kJ-N,k ^ •••^ — ^N,k [A.ZO) 

If we substitute in this expression the general formula (A. 17) for YJ^^'""* and formula 
(A. 22), we obtain 

BN,k;s = ^JV,fc;s ■ (A. 26) 

Therefore we get the general expression 

w^i ci r(^ + k-s-i 



^N,k - 1. (-1) ^ r(^ + k-l^ ^^'-^ ^^-'^^ 



We must now check the (A. 27) satisfies properties 3 and 4. Property 4 follows 
immediately: 

Yai...ak;Pi...Pk rpPi-.-Pk }_ /^o^i---o^k;fh---Pk rpfh---Pk rpai...af^ /a r)o\ 

^N,k ■^N,k — V(fc;0) ■^N,k — ■^N,k 1 (^A.ZO; 

where in the first step we used the tracelessness of Tjv,fc and in the second its 
symmetry. In order to prove property 3 let us introduce 

^r.:;r(u) ^ {y'rp^;^^^ (A.29) 

Lfc/2J 

f^^V-^u) ^ W E ^iv,fc;.^(i:;r(u) (A.30) 

s=0 



where u is an arbitrary vector. We note that 
d_ _d_ 



due. due,, e^^^) ^ ^ " ^^^^-^'^ ^ ^ ^ 
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and therefore we can write 

If) r) ^ 

T-ai •••"A: ;/?!•••& /A 

Also, from (A. 29), we obtain 

P^;^" (u) = y' P^:^^ (u) + {N + 2k-2s-2) P--- (u) ; (A.33) 
from this and (A. 15) it foUows that 

6a,aA:r'{^) = (A.34) 

and therefore property 3 is satisfied. 

FinaUy, using (2.25), we can prove the trace formula (2.16). Indeed from the 
orthogonality relations (2.10), summing over all indices we have 

^ jd^{(T) YN,k{(T)-YN,k{(T) (A.35) 

The scalar product in the r.h.s. is rotationally invariant and as such it does not 
depend on a. Choosing <t = w = (1, 0, . . . , 0) and using (2.25) we get 

^ Y^^ki^) ■ FAr,fc(w) = Y^fiw) fiN,k = ^N,k (A.36) 

A. 4 Expansions in Terms of Hyperspherical Harmonics 

We want now to discuss the convergence of the expansion (2.18). We will begin 
by showing the following result: given a generic (real) tensor T^^^^""*, the hyper- 
spherical harmonics satisfy the inequality 

[TN,k ■ YN,k{a)? < {TN,k ■ TN,k) MN,k . (A.37) 
Indeed, using Schwarz's inequality and (2.25), we get 

\TN,k ■YN,k{(J')\^ < (rjv^fc • Tjv^fc) [ Fjv^fc (cr) • Fjv^fc (cr) ] 
= {TN,k ■ TN,k) [ YN,k (w) • YN,k (w) ] 
= {TN,k ■ TN,k) J^N,k (A. 38) 

Moreover, equality in (A.37) is possible only for those a for which 

TnT' = 7 Y;^:r^ (a) (A.39) 

for some constant 7. This requires first of all T^^^^""* to be symmetric and traceless. 
The constant 7 is easily obtained squaring the previous relation: 

7^ = ^^""'i; ^^'^^ . (A.40) 



55 



Now let us consider the special case T^^^^""* = Y^^j^""'' (w) with w = (1, 0, . . . , 0) 
and k > 1. Equality in (A. 37) is possible only if 

YM:r'i^) = ir-r^w) . (A.41) 

We will now prove that if iV > 3, this implies cr = ±w. Let us firstly notice that if 
a satisfies (A.41), then every a' = Rcr with R G SO{N) such that Rw = w also 
satisfies (A.41). Now if <t 7^ w there exists an index a ^ 1 such that a" ^ 0. If 
iV > 3, we can consider rotations R in the (a,/3)-plane (with fi ^ 1) and generate 
solutions a' with cr'" assuming any value between —a" and a" . This means that 
(A.41) with «! = ... = «/; = « has an infinite number of solutions, which is 
impossible as this is a polynomial equation in a" . Thus for iV > 3 we must have 
(T = ibw. This result can easily be rephrased in terms of Gegenbauer polynomials: 
since C^^"^ ^(1) > 0, it implies that 

\C^''-\-)\ < C^''-\l) (A.42) 

for — 1 < x < 1 and k>l. 

For N = 2 the previous result is not true. Indeed in this case every cr = 
(cos vrj /A;, sin vrj /A;) with j = 1, 2, . . . , A; is a solution of (A.41). To show this let us 
notice that dim £^2,fc = 2 for all k so that there are only two independent equations 
to satisfy. Using complex indices it = l±i2 and noticing that, since S'^'^ = 8 = 0, 
we obtain 

Y^,k-^ = l^2,k<y+...a+ = ii2,ke''' (A.43a) 
Y2~k"-~ = fi2,k (T~ . . . (T~ = fi2,k e-'^^ (A.43b) 

where cr = (cos ^, sin ^). Therefore, equations (A.41) are equivalent to 

e'''^ = ±1 (A.44a) 
^-,ke ^ (A.44b) 

which proves the result. 

To discuss the convergence of the series (2.18) let us firstly notice that 

00 00 „ 

fc=o fc=o 

(A.45) 

Using the completeness relation (2.20) we get 

E /r /r = / d^{<^) \f{<^)\" , (A.46) 

which is the Plancherel identity for harmonic analysis in S^~^ . Now let us consider, 
instead of /, the function where C is the Laplace-Beltrami operator. In this 
case y^i---"* is replaced by A^^/^^'""*, and thus we obtain 

E ^N,k Ik'-"' Ik'-"' = / d^icT) \C-f{cTt . (A.47) 

fc=0 
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If now / is (7°° function, the integral is finite for all n. Thus the sum on the 
l.h.s. is converging for all n. As Ajv,fc ~ for A; — )• oo we get that, for every n, 
j^2n j:ai---(y.k j:ai---a.k _^ g ^ _j. qq_ This implies that all coefficients decrease 
faster than any inverse power of k. To prove the convergence of the series (2.18) it 
is then enough to notice that |i^jvfc — {^N,k)^^'^ and that A/jv,fc behaves for large 
k as k^-\ 

In general we can write 



J{(^ = JN;k,h ^N,k 



k,h=0 

Invariance under rotations gives 



N,k 



{R)T\ 



•fh---Ph]Si---Sh 

N,h 



I -"-J JN;k,h 



'ai...ai,;f3i...f3fi 
N:k,h 



for every rotation R G SO(N). Then, by Schur's lemma. 



JN:k,h 



ai...ak',Pi...Ph 



F, 



N,k 



SO that 



f {(T ■t) = FN,k YN,k {(t) ■ YN,k (t) 



(A.48) 



(A.49) 



(A. 50) 



(A.51) 



k=0 



Let us now discuss the properties of the coefficients -Fjv,fc in (2.27). The second 
property follows immediately from the previous discussion. We want now to prove 
that, if f (t) is positive for t G [—1,1], then |-Fjv,fc| < -Fjv,o for k > 1. Indeed from 
the definition and (A. 42) we get 



IF, 



N,k\ 



< 



< 



dt (1 - t2)(^-3)/2 f{t) 



c 



N/2-1 



Q,N/2-l 



(1) 



dt (1 - t2)(^-3)/2 /(t) 



Fn,o ■ 



(A.52) 



A. 5 Clebsch-Gordan Coefficients 

Let us now discuss the computation of the Clebsch-Gordan coefficients (2.36). 
For arbitrary completely symmetric and traceless tensors Tjv,fc, Un^i and Vw,™, we 
want to compute 



Jdn{a) [TN,k-YN,ki(T)] [UN,l-YN,li(^)] [VN,m-YN,„^i(T)] . 

Using the definition (2.3) of the hyperspherical harmonics this reduces to 
I^N,k I^N,i I^N,m I dQ{a) [a"'... a^' . . . a^' . . a""" 



rpai...ak ttPi-.-Pi T^7i...7n 
^ ■^N,k '^N,l ^N,m 



(A.53) 



(A. 54) 
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From this we see that the integral vanishes if A; + / + m is odd. On the other hand 
if A; + / + m is even we can use (A. 8). We must now compute how many scalars we 
can construct with the three tensors. It is easy to see that there is only one possible 
scalar, with the following structure: i indices of Tjv,fc are contracted with i indices 
of Vw,™, j indices of Vw,™ are contracted with j indices of Un^i and h indices of Un^i 
are contracted with h indices of Tjv,fc. Here i = {k -\- m — l)/2^ j = (m — k -\- l)/2 
and h = {I -\- k — m)/2. Of course z, j and h must be positive and this is equivalent 
to |/ — A;| < m < / + A;. We must then compute the combinatorial factor, i.e. in how 
many ways this scalar can be constructed. We find 

(k\ (l\ (m\ ,, ., k\ l\m\ 

Thus the integral becomes 

IJ'N,k f^N,l f^N,m k\ l\ ml rpai...a,bi...bh jjbi...bhCi...Cj -j^ai ...a,ci ...Cj / . 

(k+jy. ^j\h\ ^^'^ • ^^-^^^ 

Formula (2.36) immediately follows. 

Finally we want to discuss the computation of C^. Using (2.36) and (2.34) 

we get 

^N;k,l,ra " / ""(."J ^N;k,l,ra ^ N,k \" ) ^ N,l \" ) ^ N,m \" ) 

_ IJ'N,k Atjv.f Atjv.m k\ V. ml 

X / dn{a) Y^^-'^^-^^ (a) F^;/--'^---^ (^)F;--^-----(^) (A.57) 

where i = (A; + m — /)/2, j = (m — k -\- l)/2 and h = {I -\- k — m)/2. The quantity 
which remains inside the integral is a scalar; as such it is <T-independent and thus 
we can drop the integration. 

To compute the remaining contraction let us use the general expression for the 
hyperspherical harmonics given in (A. 17). Then a straightforward combinatorial 
exercise gives 

pai...a,/3i.../3fe / X ^rPl...Phll■■■lJ ( _\ t^7i •••7j «i / _\ _ 



{k;s) 

h\ I i 



(A.58) 

Note that this gives zero if s > /i or s > z. Now let us define as 

a" Y^^r''-' ia) = X, Yi;],-fr (A.59) 

It follows that 

li^m = a"^ ...a"^Y^lf"'={a) = ^ (A.60) 



m = l 



fJ'N,k 
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and therefore we obtain 

X, = (A.61) 

Thus (A. 58) becomes 

, fh\ft\ A/Vj A/jv 
■5! ■J^N,3+s 77 77 (A. 62) 

\S J \S J flN,l fJ'N,m J^N,l-h+s J^N,m-i+s 

and we get the final resuh 

2 _ f^Nk ^N,l J^N,m k\ l\ ml 

f (-i)- r(f + ^-3-i) m/A 

where M = min( [A;/2j , z, /i) with z = (k -\- m — /)/2, j = (m — A; + /)/2 and 
h = {I -\- k — m)/2. We remind the reader that (A. 63) holds only when A; + / + m is 
even and |/ — A;| < m < / + A;; in all other cases, C^. /; ; „ = 0. 

There is another way of computing C^. ^ ; Using (2.35) we can write 

CN;k,i,m = j dQia)dQ{T) [¥^4^) ■ YnA-t)] [YnA(^) ■ YnAt)] [YnA^) ■ YnA-t)] 

(A. 64) 

The integrand is only a function of <t • r. Thus, using the rotational invariance 
of the measure, we can fix one of the two spins to an arbitrary value. Let us set 
T = w = (1, 0, . . . , 0). We obtain, after integrating in c/0(t), 

Cn; k,i,m = j dn{(T) [Yj^^a) ■ Fiv,fc(w)] [Yj^^a) ■ Yj^^w)] [Yj^Aa] ■ Yj^^w)] 

(A. 65) 

and by using (2.23) we end up with 

Sn 1 

^N- klra = ^N,k J^N,l J^N,m X 

' ' ' ON 

/•I r^/^"V/^ r^/^"V/^ 

X dt{l- fr-^y^ %7TT^ %JTT^ Nn i (^-66) 
•^-1 C-fc^ (1) c-/^ (1) C^' {I) 

If one of the three indices k^l^m is fixed to some specific value this integral is 
easily done using the recursion relations of the Gegenbauer polynomials and their 
orthogonality properties [31]. In this way we have checked the general formula 
(A. 63) for A; = 1,2 and /, m arbitrary. 

A. 6 6— j Symbols 

In this section we will discuss briefiy the 6-j symbols. In dimension d = 1 they 
appear in the computation of the four-point function (not treated in this paper). 
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Figure 11: Graph showing the spin assignments in the 6-j symbol 7^jv(o, &, c; c/, e, /). 
Each vertex denotes a Clebsch-Gordan coefficient. 



and in higher dimensions they play a crucial role in deriving high-temperature 
expansions even for the two-point function. 

The 6-j symbols (also called Racah symbols) are 0(iV)-scalars defined by 



-7? ( r, h r-- rl o f \ — pcti ...a a\fh ---P d\ll c pcti ...a a\Sl ...Si,;t 
'VJVl^a, 0, C, a, e, /J — L-jV;a,(i,c '^N;a,b,e 

-,l3i...l3d;ti...tevni---'nf pSi---S 



different conventional notation for iV = 3 is 



See Figure 11 for a graphical representation. The tetrahedral symmetry which is 
enjoyed by the 6-j symbols for iV = 3 [35] is trivially true also for generic N . A 

a 6 c 1 
J d e]- 

We have not yet been able to compute a general formula for the 6-j symbols, 
but we have computed a very large class of special cases: among others, those in 
which one of the spins (say, a) takes the value 1 or 2, while the other five spins take 
arbitrary values. This class of special cases is sufficient for computing the high- 
temperature expansion of the S^~^ cr-model in general dimension d up to rather 
high order [28]. 

We begin by deriving a completeness relation for the 6-j symbols. To do this, 
let us first prove two properties of the Clebsch-Gordan coefficients. Using their 
definition in terms of hyperspherical harmonics and the completeness relation (2.20) 
we can easily prove the crossing relation 

E^ai...ap;/3i.../3j:;7i...7i pai...ap;Si...Sm;ti---tn ^ai ...ap;/3i ;5i •••5m ^ai •••ap;7i •••7; ;ei •••£« 

p p 

(A.68) 
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The second relation we need, which foUows immediately from Schur's lemma, is 

^ai...ak;Pi...Pi;ii...lm r'ai...ak;Pi...Pi;Si...S„ r ^ 7-71 •••7m ;5i •••5m /-'2 f\ pQ\ 

^N;k,l,m ^N;k,l,n ~ '^ni,n ,r J-N,m ^N;k,l,m ■ [-n-.W) 

Inserting these two relations in (A. 67) we get 

1 

Next let us compute 

1^ pSi---Spai...af,;Si...Spfji...fji;-ii...-irn pa.i...af^;fil...fiirfl..."frn A 71 ^ 

r^N,p;k,l,m — '^N;k+p,l+p,m '^N;k,l,m ■ l^-'-'-J 

Using (2.36) we get 

IJ'N,k fiN,! IJ'N,m kill ml 



Y,T^{a,'b,c-d,eJ) = j^C^.^^^^jC^.j^^^j . (A. 70) 



N,p;k,l,m 



-,Si ...Spai...a,f 

^N;k-\-p,l-\-p,m 



^(5i. .. SpQfi. .. a,/3i. .. /3ft ;(5i...(5p/3i.../3ft7i...7j;7i...7j / . 
'-'N:k+pJ+p,m ■ (^A.rZj 



where i = (m + A; — /)/2, j = (m — k)/2 and h = {I -\-k — m)/2. Now, using again 
(2.36) we also have 

^2 _ /ijv,fc+p Atjv,f+p /ijv.m (A; + p)\ {I + p)\ ml 

^N;k+p,l+p,m ~ ■ -• ..•./, x. X 



IJ.%^k+3+p(k+J + py- tljl(h + p)\ 

-,Si...Spai...a,f 

-^N;k-\-p,l-\-p,m 



^(5i. .. SpQfi. .. a,/3i. .. /3ft ;(5i...(5p/3i.../3ft7i...7j;7i...7j / . „„x 
'-'N:k+pJ+p,m ■ [-t^-IO) 



Comparing, we get 



{k+j +p)l 

'^N,p;k,l,m — 2 /; i 'M ^ 

I^N,k+P l^'N,l+p I^N,k+] + " 



k\n{h + p). 



^N:k+pJ+p,m ■ (A. 74) 



(A; + p)! (/ + p)! hi ^'>'+P''+P 

Using this result we can now compute 7^jv(p, k -\- p^l -\- p;l^k^ m). Indeed, using 
(2.36) we get immediately 

7^jv(p, k + p,l + p;l,k,m) = ^^-p ^^'^ lCN,p;k,i,m ■ (A. 75) 

I^N,p+l I^N,p+k 

Let us now derive a few other particular cases which are relevant for high- 
temperature expansions. 

Let us first consider the case in which one of the spins appearing in the 6-j 
symbols (say, a) is 1. In this case, using the tetrahedral symmetry, one can see 
that all non- vanishing symbols can be rewritten as TZ^il^k + 1,/ + 1;/, A;,m) or 
7^jv(l, A; + 1, / — 1; /, A;, m), with A;, / arbitrary and |A; — /|<m<A;-|-/in the first 
case, max(|A; — /|,|A; — / -|-2|) < m < A;-|-/in the second one. The first quantity 
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is a particular case of (A. 75), while the second one can be computed using the 
completeness relation (A. 70). Indeed we have 



7^Jv(l, A; + 1,/ - 1; /, A;,m) = -7^Jv(l, A; + 1, / + 1; /, A;, m) + 



AT, 



N,k 



2 n2 

k,k+V^N;k,l,m ■ 



(A.76) 

Let us next consider the case in which one of the spins (say, a) is 2. Using the 
tetrahedral symmetry one can rewrite all the non- vanishing 6-j symbols in one of 
the following forms: 



^N]k,Wm = 7^^^(2,A; + 2,/ + 2;/,A;,m) 
A^mI,^ = 7^^^(2,A;-2,/ + 2;/,A;,m) 



N:k,Lm 



(A.77) 
(A.78) 
(A.79) 
(A.80) 



with k and / arbitrary, m < / + A;, and m > |/ — A;| for A^j^\i^ and A^^\i^^ 

m > max(|/ — A;|, |/ — A; + 2|) for A^^],.^i^^, and m > max(|/ — A;|, |/ — A; + 4|) for A^^],.^i^^. 
Using the completeness relation (A. 70) we can rewrite the last two quantities in 
terms of the others. Indeed 



N:k,Lm 



1 



N:k,Lm 



/l(2) J n2 

■^N;k,l,m^ »/■ '^N;2,l,l+2'^N;k,m,l 



N:k,l,m 



k,rri 



A%\ira ^ particular case of (A. 75). To compute ^jv-fc/m' use (2.36) to 

get 



+ ^2 

A/jv,/ 
1 

(2) 



(2) a(2) j ^_r2 n2 

■^N;l,k-2,m »r 



n2 n 
^N:2,Ul^- 



N:k,m,l 



(A.81) 
(A.82) 



(2) _ 2A; /^Af,2 l^'N,k l^'N,l ^af3-n...-ti;a5i...5k-i;vi-Vmr'^fi--n;f3Si...5k-i;vi-Vr, 



^N:kd,m "7,1 2 ^Nd+2,k,m 
* + A«JV,fc+l I^N,l+2 

Then, using again (2.36) and (A.22)/(A.23)/(A.26) we get 

a(2) _ 2k /i^^2 f^N,k f^N,l ^af3-n...-fi;a5i...5k-i;vi' 

■^N:kd,m — - - o C-A 



k,rri 



(A.83) 



k + l I^N,k+i I^N,l+2 



-'N;l+2,k,m 



fJ'N,k fJ'N,l ^f3ji...ji;5i...5k_i;rii...rim . 

■ '--JV;/+l,fc-l,m "T 



I^N,k-l I^N,l+l 



/J, -r \ 4 _J^^Njkcf5Si pii---ii\S2---Sk-ivni---'n-n 

[K — l)/iN^k;l C '-^Nd,k-2,m 



fJ'N,k-2 



(A.84) 



[where Ajv,fc;i is defined in (A. 16)], and thus 



A 



(2) 

N:kd,m 



N{N + 2) 



k 



I^N,k I^N,l 



k + I fJ'N,k+l fJ'N,k-l fJ'N,l+2 



l+l,k-l,m 



fJ'N,k-2 
I^N,k-l 



JV,2;/,fc-2,' 



(A.85) 
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B Finite- Size- Scaling Functions for the Univer- 
sality Classes (4.84) 

In this appendix we want to study the finite-size-scahng functions for the one- 
parameter family of universahty classes (4.84): this family is parametrized by a 
real number B G [0, +00], and interpolates between the iV-vector universality class 
{B = 0) and the RP^~^ universality class {B = 00). In particular, we want to study 
the asymptotic behavior in the perturbative regime (7 — )• 0), and show that in the 
even-spin sectors (k = 2,4,...) the finite-size-scaling functions are independent of 
B modulo nonperturbative corrections of order roughly e~'^ (We will succeed 
here in doing this only for A; = 2, but we conjecture that it is true for all even k.) 

The basic idea can be seen in the simple case of the partition-function scaling 
function ^^^^(7), defined in (4.77). We have 



zPir^B) = J2^N,ie-™ (B.la) 

00 00 

= E e-"'"^-' + e-"''' E e-"'"^-' (B.lb) 



/ = / = 

/ even / odd 



1 _L g-7-B _ I _ g-7-B _ 

—l^ZUl) + -^^Z^ij), (B.lc) 



where we have defined 



Z±(7) ^ E (±1)' ^N,i e--'-^' . (B.2) 

It is easy to see that Z^{j) is of order 7-(^-i)/2 as 7 — )• 0: roughly speaking, for 
small 7 the sum over / can be replaced by an integral. On the other hand, we 
shall show that Z^(j) is exponentially small as 7 — )• 0: more precisely, it is of 
order e~'^^/'*"'7~(^~2). Thus, the i?-dependence of Z^]^\j;B) is given by the trivial 
prefactor (1 -|- e~"'^)/2, up to nonperturbative corrections of order roughly e~'^ 

A similar result will be shown for the numerator of the susceptibility scaling 
function xjv fc(7j ^) [^^^ (4-79)] for k = 2, from which it will follow that xjv^2(7j ^) 
is independent of B modulo nonperturbative corrections of order roughly e"'^^ 

The crux of the matter will thus be to control the behavior of ^^(7) [and the 
analogous numerator functions] as 7 — )• 0. For N = 2 this is a simple consequence 
of the Poisson summation formula [(B.6) below]: 

^^=2(7) ^ E = - E e-^^^/"^"^ (B.3) 

/= — 00 \^/ m = — 00 

00 / — \^/^ 00 -I 

^^=2(7) - E (-1)'^-^' = - E e-(-^/-)(-+i)^ (B.4) 

/=-oo \7/ 



-00 



63 



Likewise, for = 3 the Poisson summation formula controls Z (though not 

/= — oo \^ / m = — oo 

(B.5) 

We shall derive the analogous identities for general iV, by two alternative methods: 

(a) by developing generalizations of the Poisson summation formula (Sections B.l 
and B.4); and 

(b) by differentiating the identities for N = 2 or N = 3 with respect to 7 [e.g. 
(B.42)]. 

B.l Generalized Poisson Summation Formulae 

Let / G S(R) [i.e. / is a function of one real variable that is infinitely differen- 
tiable and that together with all its derivatives vanishes at infinity faster than any 
inverse power], and define f(t) = Jf^ e~^^^ f{x) dx. Then we have the well-known 
Poisson summation formula 

00 00 

E m = E /(2vr0 (B.6) 

k= — oo /= — 00 

and its (less well-known) one-sided generalization 

00 00 ■ °? , 

/(O) + 2 E f(k) = E /(2vr0 + cot - f{t) dt , (B.7) 



where P denotes Cauchy principal value at each of the singularities of the inte- 
grand. For a proof, see [36, pp. 31-32 and 64-65]. Here we will prove the following 
generalization of the one-sided Poisson formula: for any real iV, 

1 7 ^ 

Y.^N,kf{k) = — J KN{t)mdt , (B.8) 

A;=0 — 



where 



1 _|_ g«(*+«) 

KN{t) ^ limj^—-^^^-^ (B.9a) 

^ ^-m-2)/2^.^ (B.9b) 

40 (-2zsin^)^-i ^ ^ 

= e-''^^'-'^' LN{t) (B.9c) 
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is a well-defined distribution in iS'(R). Note also the recursion formula 



1 



N{N - 1) 



dt 



2 



LN{t) . (B.IO) 



For N integer > 2, we can make (B.8)/(B.9) more explicit: 

oo JV-2 oo 1 1 _|_ p«* 

j:^N,kf(k) = E ^N,n E /^"^(SvrO + ^P/ _\,N-i (B.lla) 



N-2 oo m 

ra = /=-oo 
, + even 



t=27r/ 



— oo ^ 

where the cjv,™ and cjv,™ are Laurent coefficients (see below), and of course p'"''^ 
denotes the n^^ derivative of /. For N = 2 we will have C2,o = C2,o = 1, so that in 
this case (B.ll) reduces to (B.7). 

We start from the well-known identity 

^ T(N + k) z'' , ^ ^ , , 

valid for complex z in the disc < 1. (If iV is noninteger, we of course take the 
branch that equals 1 at 2; = 0.) Using the expression (A. lb) for A/jv,fc, it follows 
immediately that 

00 

Fn{z) = J2^N,kz' = (1-Z)-^(1-Z2) 



k=0 



(1 -z)-(^-^)(l + z) . (B.13) 



We will use this identity to construct a "complex- variables" proof of (B.8)/(B.9) and 
(B.ll). [It would be interesting to know whether there is a simple "real- variables" 
proof, as indeed there is for (B.6) and (B.7): see e.g. Ref. [36], pp. 31-32 and 64-65.] 



®'^The existence of (B.9) in the sense of tempered distributions is a consequence of the following 
theorem [37, Theorem 2-10, pp. 62-63]: If / is analytic in the strip < Imz < R and satisfies 
there the bound \f(x + iy)\ < C(l + \x\P)y~'' for some C,p,r < co, then limyiofix + iy) exists 
in 5'(R). Sketch of Proof: For a test function g G <5(R), define h{y) = J f(x + iy)g{x)dx. 
We can compute the derivatives of h using the analyticity of / and integration by parts: h^"'\y) = 
(— i)" / f{x -\- iy) g'-"\x) dx. It follows that \h'-"\y)\ < Cn \\g\\p,n y~^ for a Schwartz norm || • 
Starting from n > r and using the fundamental theorem of calculus, it is easy to show that limj,|o h(y) 
exists, with uniform bounds in terms of a Schwartz norm of g. Q.E.D. See also [38, Section 12.2, 
Corollary 4, p. 192] for a similar result. 
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Let us begin by assuming that the function /, in addition to lying in iS(R), 
satisfies the bounds \f ''Hx)\ < C^e"^!^! for some constants (7„ < oo and S > 0; later 
we will relax this assumption. It follows from this that f(t) is analytic in the strip 
I Imt| < S; moreover, in this strip / vanishes faster than any inverse power of | Ret| 
when I Ret I — )• oo. So we can use the representation 

oo 

f(k) = — f e^''mdt, (B.14) 

— CO 

where the contour of integration runs slightly above the real axis (say, at Imt = e 
with < € < S). It follows that 

oo -1 oo ^ 

J2MN,kf(k) = — J2^N,k / (eyf(t)dt. (B.15) 

k=0 k=o _4, 

This joint sum/integral is absolutely convergent (since |e'*| = < 1 and f(t) 
decays rapidly at infinity), so we can interchange the summation and integration. 
Using (B.13), we obtain 

t^N.fik) = ^J ];%_^ mdt , (B.16) 



fc=0 



where the integration still runs at Imt = e. Since the value of the integral is 
independent of e (for < e < (5), we can trivially take e J, 0, thus proving (B.8)/(B.9) 
for functions / satisfying the above restrictions. 

It is easy to remove the assumption that / and its derivatives decay exponen- 
tially. Just apply the foregoing result to fa{x) = f{x)e~"^ and let a J, 0. Then /„ 
equals / convoluted with a Gaussian (47ra)-i/2e-*'/4«, and this G aussian tends (in 
the sense of distributions) to a delta function as a J, 0; therefore, /„—)•/ in S(R) as 
a J, 0. In particular, the right-hand side of (B.8), taken on /„, tends as a J, to its 
value taken on /: this is an immediate consequence of the fact that A'jv(0 defines 
a distribution in iS'(R). On the other hand, the left-hand side of (B.8) converges to 
its a = value by virtue of the dominated convergence theorem. 

Let us now assume that N is an integer > 2, and let us again temporarily 
assume that / and its derivatives decay exponentially. Then the integral (B.16) at 
Imt = e can be written as the half-sum of the integrals taken over Imt = ibe plus 
the half- difference. Now the half-sum is, by definition, precisely the principal- value 
integral in (B.lla)^^; on the other hand, the half- difference is —iri times the sum of 
the residues at the poles t = 2ttI {I integer). Using the Laurent expansion 



1 + e'* 



V 

(JV-1) 



(l_e»-t)iv-i - ^ «JV,m^" (B.17) 



®^Our "principal-value integral" is the same as the "canonical regularization" of Gel'fand and 
Shilov [36, Sections 1.3 and 1.4]. Note, in particular, equations (6), (7), (12) and (13) in Section 1.4.4 
(pp. 94-95). 



66 



around t = (and of course an identical expansion at each pole t = 27r/, / integer), 
we obtain (B.lla) with 

CN,n = -- j • (B.18) 

We can see that the leading term in the Laurent expansion is ajv,-(Af-i) = 2z^~^, 
and hence the highest-derivative coefficient in (B.lla) is cjv,Af-2 = —i^/{N — 2)!. In 
particular, for N = 2 we have C2,o = 1, so that (B.lla) reproduces (B.7). 
Equivalently, we can use the Laurent expansion 



2cos^ 



2 



(-2z sin 1)^-1 



2/ m = -(JV-l) 

around t = (and of course an identical expansion multiplied by ( — 1)^' at each 
pole t = 27r/, / integer). Note that ajv,m 7^ only when m + iV is odd. We therefore 
obtain (B.llb) with 

CN,n = -- j • (B.20) 

z nl 

Note that cjv,™ 7^ only when n+iV is even. The leading terms are ajv,-(Af-i) = 2i^~^ 
and hence cjv,Af-2 = —i^/{N — 2)\. From (B.IO) we can derive the recursion relation 



1 



N{N - 1) 



(m + 2)(m + l)ajv,„+2 + ( — - — j o^N. 



■iV - 2\2 ^ 

(B.21) 



which together with the initial conditions a2,_i = 2i and 0,3^^2 = —2 yields all the 
coefficients. Unfortunately, we have been unable to find a closed-form solution for 
this recursion relation. 

The assumption that / and its derivatives decay exponentially can be removed 
as before, using the fact that both terms on the right-hand side of (B.lla)/(B.llb) 
define distributions in iS'(R). 

If iV is a real number < 2 (not necessarily integer), we can rewrite the kernel 
K]y(t) in a somewhat more explicit form. Note that Re(l — e'^*"*"")) > for all t, 
hence | arg(l - e'(*+"))| < 7r/2. It follows that 



cos * 



I 2 1 

where 

N -2 

Lp{t) = — X (t mod 27r) (B.23) 

and t mod 27r is taken to lie in the interval [0,27r). Since iV < 2, (B.22) defines 
a locally absolutely integrable function, hence is unambiguous as a distribution. 
Equivalently, we can write 

r — i^^t T,^ /,\ o2-iV i(JV-l)7r/2 | COS ^ | -^,(4) {T>nA \ 

LNlt) = e 2 Ajv(t) = 2 ' ' . ^ e ^ , (B.24) 

I 2 I 
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where 

N - 2 

ib{t) = if{t) + t = 7r(iV - 2) [t/27rj (B.25) 

and [xj denotes the largest integer < x. 

For N > 2 these formulae are ill-defined because A'jv(0 nonintegrahle sin- 
gularities at t = 2ttI {I integer). For noninteger iV > 2, explicit formulae can be 
obtained by using the recursion formula (B.IO) [starting from (B.24) at some N <2] 
together with integration by parts. (More precisely, integration by parts is how one 
defines the derivative of a distribution!) 

Finally, let us go back to (B.9b,c) and note an interesting property of the kernel 
-Ljv (valid for all N): we claim that if we decompose L]\j{t) into its symmetric and 
antisymmetric parts around t = vr, 

L%{t) = \[LN{t)±LN{2T^-t)] , (B.26) 

then the symmetric part L^{t) vanishes on the interval < t < 27r (i.e. it is 
supported outside this interval). PROOF: The numerator cos is obviously anti- 
symmetric around t = vr in the limit e J, 0. As for the denominator, the function 
values sin and sin (^^~^*)+'^ belong to the same Riemann sheet of the function 
-;JV-i provided that < t < 27r (and not otherwise), so that in this case they tend 
as e J, to the same point on the Riemann surface. Therefore, the denominator is 
symmetric around t = vr in the limit e J, 0, for < t < 27r (and only there). Q.E.D. 

Of course, the same argument can be made around any point t = {21 -\- l)7r, / 
integer: the symmetric part vanishes on the interval 2ttI < t < 2tt{1 -\- 1). 

This symmetry /support property is of particular relevance in case the function 
e~' 2 * f{t) is symmetric around t = tt (as will be the case in our application below). 



B.2 Some Generalized Theta Functions 

Now we want to apply the generalized Poisson summation formulae to analyze 
the asymptotic behavior as 7 — )• of some generalized theta functions. Let us define 

00 

Z^^eAl) = J2^N,ke'''e-^^'+"^' , (B.27) 

k=0 

which of course is periodic in 6 with period 27r. Applying (B.16) with f{x) = 
g.feg-7(^+«)2 obtain 



1/2 



00 



1 / TT \ -,0 f 2 COS I 



27r \ j J " J (-2zsin|)^- 



-e^i»-^)^e-^'-'^"/^''dt (B.28b) 



-00 



1 fTrV^' ... 2z 7 1 d 



2Tr \-fJ ^ iV-2 y (-2zsin|)^-2 (/t ' ' 

(B.28c) 
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where the contour of integration runs at Imt = e > 0; here (B.28c) is obtained from 
(B.28b) by integration by parts. Note that the formula becomes shghtly simpler in 
the case a = (N - 2)/2. 

Let us consider first the case of N integer > 2. As before, the integral (B.28b) 
can be written as a principal- value integral plus —iri times a sum of residues. To 
compute the residue contribution, we use the Laurent expansion (B.19), yielding 



residue contribution to Zjv,e,a(7) 

JV-2 



^\l/2 JV-2 oo in 



7 / n ' , dt 



t = 2TTl 

(B.29) 



The sum over / is absolutely convergent, uniformly on compact subsets of the half- 
plane Re 7 > 0, thanks to the rapid decay of e~(*~^) as t — )• iboo. Moreover, as 
7 — )• this sum is dominated by its leading term(s), namely the one(s) for which 
|27r/ — ^1 is smallest. 

Concerning the principal- value integral, we first note that in certain cases it 
vanishes by symmetry: If N is an integer, then cos |/(sin |)^~"^ has parity ( — 1)^""^; 
it follows that the combination 



ZN,e,ail) + (-ire-(^-2-«)^Zjv,_,,jv_2-„(7) (B.30) 



is given exactly by the sum of residues. In particular, in two cases Zjv,e,a(7) itself 
is given by the sum of residues (B.29): 



(a) a = {N - 2)/2, 6 = 0, N 



even: 



' \l/2 Af-2 oo vn 



\ / / n=0 /= — oo 



As 7 — )• 0, this sum equals its / = term (which is of order 7 ^^Z^) up to 
nonperturbative corrections of order e~'^'^ For N = 2 this reduces 
to (B.3). 

(b) a = (N — 2)/2, 6 = TT, N integer (even or odd): 

/^\l/2 JV-2 00 in 

\ ' / n=0 /= — 00 



n=0 l- 



-00 



t=27r/ 

(B.32) 

As 7 — > 0, this sum is exponentially small, of order e~^'^ For N = 2 
(resp. iV = 3) this reduces to (B.4) [resp. (B.5)]. 
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If one makes the change of variables t = s + tt and then uses the oddness of the function 

2/(cos f ; 



sin |/(cos ^ (for all integers N), one obtains the same combination (B.30) 
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Next let us consider the general case of N real (not necessarily integer). Sym- 
metrizing (B.28b) around t = vr, we find 



W f^A _L JttN -i{N-2-a)e rpr f \ 

1 /vr 



,1/2 oo 



27r \ 7 / 

— oo 

(B.33) 

where L^{t) is the symmetric part oi L]\f{t) around t = vr [cf. (B.26)]. As discussed in 
the preceding subsection, L^{t) is supported away from the interval < t < 27r, and 
is the p^^ derivative of a polynomially bounded function, where p = max([iVj,0). 
The following lemma then implies that (B.33) and all its derivatives with respect to 
7 are exponentially small whenever 9 mod 27r: more precisely, the n^^ derivative 

of (B.33) is bounded by const X ^~(?'+ 2 "'"^"^e"®^ where 

= min|^ - 27rk\ . (B.34) 
fceZ 



In particular, if a = (iV - 2)/2 and < ^ < 27r, then Re[e'(~'^~t^) (^'^] 
IS as 7 J, (7 real). For 9 = ir, this says that Zat _ n-2 (7) is 

0(7-^6-^'/^"'). 

Lemma. Let L(t) be a tempered distribution on R, supported on |t| > 0, and 
define 

00 

F{j) = 7-^/2 J L{t)e-'"l*^ dt (B.35) 

—00 

for Re 7 > 0. Then there exist constants p and (7„ such that 

|F(")(7)| < (7^|^|-(?'+|+2-)e-e2Re(i/47) (B_3g^ 
for n > and (say) I7I < 1. 

Proof. For some p > 0, L is the p^^ derivative of a polynomially bounded function 
/, i-e. \f{t)\ < C{1 + |t|"). Thus, for each n > 0, 

FW(7) = (-1)^ Jm ^ 1^ (7-^/^ e-V4.) . (B.37) 
—00 

Clearly 



^-l/2g-tV47 



atf ^7 

Integrating over t then proves the lemma. 



< const(n) X (1 + |t|f+2-) |^|-(f+2+2-) e"*' ^^^(i/^"') . 

(B.38) 
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Finally, let us consider the cases in which the principal- value integral (for N 
integer) or the integral over < t < 27r (for N generic) does not vanish. In these 
cases we can obtain an asymptotic expansion for Zjv,e,a(7) in powers of 7, by the 
usual method of expanding the integrand of (B.28b) around its peak at t = ^. For 
6^0 (mod 27r) the result is as follows: 



(7) 



-i{N-2)e/2 gi7r(JV-l)/2 2"^-^ 



COS 



(sin 1)^^-1 



(3N -4 

1 + 



N 



a 



I a 



N 



N{N - 1) 



cot' 



+ 0(7^) 



6 6 

tan - + (iV - l)cot - 



(B.39) 



This expansion can be proven rigorously by cutting Ljv (using a smooth partition 
of unity) into a part supported on the interval e < t < 2tt — e (here we suppose 
< t < 2tt) and a part supported on the union of intervals t < 2e and t > 2tt — 2e. 
The integral over the first interval is then an ordinary integral of smooth functions, 
and the asymptotic expansion can be controlled by standard techniques; while the 
integral over the second region is exponentially small by virtue of the Lemma above. 

An alternative way of deriving these formulae is to use a recursion formula 
yielding Zjv+2 in terms of Zjv; in this way, all integer values of iV (> 2) can be 
handled by differentiating the cases N = 2 and iV = 3 with respect to 7, while all 
noninteger values of iV (> 0) can be handled by differentiating one of the cases in 
the interval < iV < 2. The basis of this approach is the simple identity 



Af - KN + k-2) 



N{N - 1) 



(B.40) 



[see (A. 2)]. It follows from this that 



JN+2,e,a + l 



(7) 



N{N - 1) 

In particular, for a = {N — 2)/2 we get 



— + 2i[a 
0(7 



^jv+2,e,f (7) 



N{N - 1) 



d_ 

c/7 



N 



d_ 

lie 



a 



N-2 (7) 



■^jv,e,a(7) • 

(B.41) 

(B.42) 



Clearly, if Z 



1(7) is exponentially small together with all its derivatives, then 



the same holds for -^jv+2r e jy+2r-2 (7) for every positive integer 



r. 



l/2^-(t-e)V47 



Remark. The equivalence of (B.IO) and (B.42) comes from the fact that 7"^ "e 
is a solution of the heat equation, hence d/ and c/^/c/t^ act identically on it; and 
since these two operators commute, the same holds true for multiple applications 
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of these operators: 



d\ Id \ / 1/2 (t_e)2/4^\ _ / 1/2 (t-e)2/4^ 



= [^y^ (,-/^e-(-^)V^.) . (B.43) 

B.3 The Part it ion- Function Scaling Function Z^^^j^B) 
Recall from (B.1)/(B.2) that 

zS\r,B) = ^^^^^(7) + ^—^Z^il), (B.44) 



where 



Z±(7) ^ E (±1)' ^N,i e--'-^' . (B.45) 



N 



ow 



Aiv,/ ^ KN + l-2) = [1 + ^) - ■ (B.46) 

So the functions Zj^ are precisely generalized theta functions of the type considered 
in the preceding subsection; indeed, they belong to the "simple" case a = (iV — 2)/2: 

Z+(7) = e"'(^-^)^/^Z^^„^^(7) (B.47a) 

ZM = e-(^-'r/^z^^^^j^{j) (B.47b) 

It follows immediately from the results of the preceding subsection that 

Z^n\t,B) = ^-^f^ZUl) + 0{e-'^'-j-(^-ly), (B.48) 

as claimed. 

Remark. The duality formula (B.3)/(B.4) for ordinary theta functions is a 
special case of a modular transformation, and is connected with the theory of el- 
liptic functions [39] [40, Chapter 13] [41]; it also has applications in string theory 
[42] . We wonder whether the corresponding formulae for integer iV > 3 are telling 
us something deep about the Riemannian geometry of the sphere S^~^. We are 
intrigued by the fact that the generalized theta functions arising from Zj^ fall pre- 
cisely into the "simple" case a = (N — 2)/2 — it can't be a mere coincidence! And 
we wonder why there is a convergent duality formula for Z^ for all integer N [cf. 
(B.32)], but for Z^ only for even N [cf. (B.31)]. Is this perhaps related to the fact 
that — / G SO(N) for N even but not for N odd? Or to the fact that the groups 
SO(N) fall into different families of the Lie classification for N even and N odd? 
And can our results be generalized to symmetric spaces other than S^~^? 
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B.4 Some More Generalized Poisson Summation Formulae 



To handle the numerator of the susceptibihty scahng function, we wiU need 
to study sums of the form ^N,k R{k) f{k)^ where i? is a rational function and 
/ G iS(R). Unfortunately, such sums are not covered by the generalized Poisson 
formulae of Section B.l: the trouble is that R is typically not a smooth function 
on all of R, so it cannot be absorbed into /. Instead, we shall derive some further 
generalizations of the Poisson summation formula, in which R is absorbed into the 
kernel K^. 

Let, therefore, i? be a rational function of the form 

where P is a polynomial. Let ko be an nonnegative integer, chosen large enough so 
that none of the fii are equal to an integer < —ko. (That is, R does not have any 
poles at integers > k^. When Ajq = we shall omit it from the notation.) We shall 
then prove a formula of the form 

1 7 ^ 

Y,^N,kR{k)f{k) = — KN;R;ko{t)f{t)dt. (B.50) 

j~r zvr J 

In fact, the derivation is virtually identical to that in Section B.l: We introduce 
the function 

oo 

FM;R;ko(z) = J2 ^N,kR(k)z' . (B.51) 

k=ko 

This series converges in the disc \z\ < 1, but FN-R-ko then has an analytic contin- 
uation to the whole 2;-plane cut along the ray [+1, -|-oo).^"^ In particular, the only 



®^This follows from [40, Theorem 11.1.3, pp. 41-43], which states the following: Let fi he a finite 
complex measure on [0,1], with |/u({l})| < p, and let //^ = fx'' d/i^x) be its moments. Let ki be any 
integer such that J x^'^ \dii{x)\ < p (such an integer always exists). Let G be an analytic function in 
the disc of radius p centered at the origin. Let /(z) = XlfcLfci z*^ be a function having an analytic 
continuation to a domain A which is starlike with respect to the origin. Then the function defined 
by XlfcLfci G{pk) Ck likewise has an analytic continuation to A. 

We shall apply this theorem as follows: Let A be the cut plane C \ [+1, -|-cxd). Let 

oo 

/(Z) = ^WiV,^^' 

fci-l 

= (l_^)-(iV-l)(l + ^) _ Y.^N.kz'' , 

k = 

where we will choose ki later. Let dp(x) = dx, so that pk = l/(k + 1). Let 

G{x) = 



l + (P-l)x 

so that G(pk) = l/(k + [3)" . To apply the theorem, it suffices to take fci > |/3 — 1| — 1. But then we 
can add in "by hand" the terms ko < k < ki — 1, provided that none of these values of k equals —[3. 
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singularity of FN-R-ko on the unit circle is at 2; = 1, and the growth of FN-R-ko this 
singularity is approached is bounded by a polynomial in |1 — 2;|~^. We can therefore 
introduce the distribution A'jV;_R;fco t>y 



KN;R-M{t) = lim FN;R-M{e'^'+''^) 
40 



(B.52) 



and its only singularities are at t = 27r/ (/ integer). The proof of (B.50) then follows 
exactly as in Section B.l. 

If A;o = 0, then Fn-r can be written explicitly in terms of the generalized hyper- 
geometric function g+i-Fg (defined e.g. in [31, formula 9.14.1, p. 1045]): 



Fn-r{^) = P[z 



dz 



n /3-'K+i^,(iV, . . . , /3i + 1, . . . , + 1; ^) 



\i=i 



+ 2)-M ,+iF,{N, /3i + 2, . . . , + 2; /3i + 3, . . . , + 3; z) 



\i=l 



(B.53) 

For 5 = 1, this special case of the 2F1 corresponds to incomplete beta functions. 

We want now to derive some general properties of these functions. Let us thus 
introduce 



T{N + k) 



v{N) k\{k + fiy 

/3-%+iF,(iV,/3,...,/3;/3 + l,...,/3 + l;z) 



(B.54) 



(Here N is fixed, so we suppress it from the notation.) By making a shift k ^ k-\-l 
in the sums, it is easy to derive the recursion relation 

UP-Z) = ^-^(l-zf-'' + (l-^^]zf,i^ + l;z) 



/3 



g-l 



{N-l)zJ2^^-^-'M^ + l;z). 



(B.55) 



n=l 



Using this formula it is easy to get Fn-r for the simplest nontrivial case R(x) 
l/(x + /3): 



°° z 1 



z{j3 + l- N) 



+ z 



/3 + 1 

2 (iV - l)(iV - 2 - 2/3) 



/3(/3 + l) 



(1-.)-^ 

Vi(/3 + 2;^). (B.56) 



This proves the claim for R(x) = l/(x + [3)" . A general denominator Q{x) then follows by expansion 
in partial fractions, and a general rational function R{x) = P(x)/Q(x) follows by application of the 
differential operator P(zd/dz). Q.E.D. 
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A simplification occurs for /3 = (N — 2)/2: the last term vanishes, and we have the 
explicit formula 

In the following we will be especially interested in the value of /g(/3; z) at z = —1. 
For q = 1 a general formula can be obtained for /3 = {N/2) + integer. The starting 
point is Kummer's formula^"* 

2Fi(a, 6; 1 + a - 6; -1) = 2-" , T{1 + a - b) ^ 

Setting a = N and b = iV/2, we get 

SO that 

/.(?;-!) = #^^- IB.COI 

Using the recursion relations (B.55), it is then possible to compute —1) for all 

13 = (N/2) + integer. 

In the following we will use two specific functions: 

oo k 

Un{z) ^ Y.^N,k-T—w^T^ (B-61) 

fc=o y'^ + 2 ) ~ ^ 

These series are well-defined provided that iV 7^ 4, 2, 0, —2, .... Simple algebraic 
manipulations yield 

Un{z) = \ [A(f -2-z)- {z' + l)A(f ; .) + .Vi(f + 2; z)] . (B.63) 
Using now the recursion relation (B.55) forward and backward to express everything 



in terms of /i(y; 2;), we obtain (specializing for simplicity to 2; = — 1) 



*See [43, p. 107, equation (47)]. See also [44, p. 50, Theorem 8.6c 
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In complete analogy we can compute Vn( — 1). Using the recursion relations (B.55) 
we get 



V( U - nN-l){N-3) ^ 



(iv-2)(iv-4) r(^) ■ 



(B.65) 



Notice that in principle one would expect here also a term proportional to /2(y; — 1); 
but for the specific combination which appears in Vn( — 1)^ the coefficient of this term 
vanishes. Notice, finally, that 

(iV-2)(iV-4)1/Ar(-l) = 2(iV - 3)t/jv(-l) . (B.66) 

A key cancellation in Section B.6 will rely on this identity but not on the specific 
values of Un( — 1) and Vn( — 1). 

B.5 Some More Generalized Theta Functions 

Let us now introduce some more generalized theta functions, which will play 
an important role in our treatment of the numerator of the susceptibility scaling 
function xjv^2(7j define, for Re 7 > 0, 

UnaM ^ Y.^N,ke^' \,^ , N^,,_, (B.67) 



k=ko \"' \ 2 

Wo(7) ^ J:^N,ke^'' (B.68) 

Here ko is an integer; if a = {N — 2)/2 is an integer, then we require that ko > 1 — a 
in order to avoid zeros of the denominators. We will thus take Ajq = except when 
N is an even integer < 4. When Ajq = we omit it from the notation. 

We remark that the functions U^^e and V^^e satisfy recursion relations identical 
to (B.42). 

Applying (B.50) with f{x) = e'^^e~"'(^+") , we obtain 

^ / \l/2 00 

UMl) = ^(H J KS'\t)e^"('-'K-('-'^'/'-dt (B.69) 

—00 

^ / \l/2 00 

VMl) = ^(H J KS'\t)e^"('-'K-('-'^'/'-dt (B.70) 

—00 

with the obvious kernels A'^^ and K^f^K For 9 (mod 27r) we can then obtain 
an asymptotic expansion of U^^ei'^) and VAf,e(7) in powers of 7, using the Lemma of 
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Section B.2 to control the contribution of the singularity. At zeroth order we have 

limUNAl) = UNie'') (B.71) 

limVNAl) = VNie'') (B.72) 

and in particular these limits are finite. For 6 = tt we have calculated these limits 
in (B.64)/(B.65). 

When 6 = TT — which we will assume henceforth — much more can be said. 
The simplest approach is to use the differential equations 

^ -Zj,^,^N^{j) + 'f:A'M,ke'''e--^'+'^^' (B.73) 

^ k=0 

(^ + l)Wo(7) = -UNAkoil) (B.74) 

to reduce the problem to known results for Z^^r _ n-2 . One can write immediately 
the solution of (B.73)/(B. 74): 

7 

UN,eM = e-^UM,e-M(0) - Je-^^-^'^Z^^,^^.JY)dl' (B.75) 



7 

VN,e;koil) = e-'VN,e;kM - j e-(^'-''^UN,e-M{i)di (B.76) 



For iV 7^ an even integer < 4, we can take Ajq = and use the fact that Zj^ ^ n-2 (7) 
is exponentially small as 7 J, 0, to get 

Un,-k{i) = ^jv( — + exponentially small terms (B.77) 
"l/jY,7r(7) = [VAr( — 1) — ?7jv( — 1) 7]e~"' + exponentially small terms (B.78) 

where ?7jv( — 1) and Vn{ — 1) have been calculated in (B.64)/(B.65). 

Next let us treat the case iV = 4, taking ko = 1. We have the initial conditions 
^4,7r;i(0) = — I and V4^,r;i(0) = ~ • Usiug the fact that Z^^t^A^) exponentially 
small, we get 

^4,7r;i(7) = (y - ^ "I" 7^ ~^ exponentially small terms (B.79) 

^4,7r;i(7) = ~ 24 4^ ~ 2^^) ^-^P^^^^^^^^^y small terms 

(B.80) 
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Finally let us treat the case iV = 2, taking ko = 2. We have the initial conditions 
^2,7r;2(0) = | and V2,7r;2(0) = " |- By the same logic we get 

^2,7r;2(7) = "'" ~ (2 ^-^P^^^^^^^^^^ Small terms (B.81) 

exponentially small terms (B.82) 

Formulae (B.79)-(B.82) can alternatively be derived from (B.77)/(B.78) and 
(B.64)/(B.65). by taking the limits iV — > 2,4 starting from noninteger N. 

B.6 The Susceptibility Scaling Function x%\{l'-,B) 

Now we want to prove a formula for the numerator of the susceptibility scaling 
function xjv fc(7i ^) — that is, for the sum appearing in (4.79) — analogous to that 
proven in Section B.3 for the partition-function scaling function Z^]^\j; B). We 
conjecture that such a formula is true for all even A;, but we shall prove it here only 
for k = 2. We define 

l,m=0 ^N;l,m 

(Note that the properties of the Clebsch-Gordan coefficients guarantee that, for k 
even, / and m in this sum have the same parity. In particular, Ajv; i,m = ^N,m — ^N,h 
independent of B.) We shall prove that X^2{l) exponentially small as 7 J, 0, so 
that the numerator in xjv^2(7i ^) becomes simply 

X^2il) + exponentially small terms . (B.84) 

From this and (B.48) it follows immediately that 

Xn\{1'iB) = TT _|_ gjjponentially small terms . (B.85) 

Nn,2 7^w(7) 

In particular, xjv^2(7j independent of B modulo exponentially small terms. 



We start by rewriting Xj^2{l) 
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^i^,2(7) = E(-l)' 

/=0 



•-JV; 2,/,/+2 '^N; 2,1,1-2 



oo 



AjV; 1^1 + 2 AjV; 1-2,1 



(B.86) 



®®In the final term we have / = m, and so we must use the comment in footnote 44 to resolve the 
ambiguity in (4.79). 
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where we set C^. 2 ; ;_2 = if / = 0, 1. We shall deal with these two sums separately. 

Let us suppose first that iV 7^ 2, 4. For the first sum in (B.86), using the formula 
(2.41) we get 



1=0 

NjN + 2) ^,(^.,2)2/4 
16 



<N - 3)^ - ^ + 2iV - 1) V^.^l) 



(B.87) 



We can then use (B.77)/(B.78) to get 

- ^^^Z^ e<(N-2fl^ _ 2)(iv - A)Vj,{-l) - 2{N - ?,)Uj,{-l) 



32 

- (iV - 2){N - 4)7t/jv(-l)] + exponentially smaU terms . (B.88) 
Using the identity (B.66), we finally get 

_ N{N-A){N ~4) g7(JV-2)V4^^^^_^^ ^ exponentially smaU terms . (B.89) 
32 

For the second sum appearing in (B.86), using (2.42) we get 

iV2 -4 



^^.{N-2fl, V(-l)'A^^. /^+ ."^ \ ^ e-^('+^)- 



7 e 



7(JV-2)V4 



d [N 



1 



Un-A.i) ■ 



(B.90) 



Using (B.77) we get 



iX^^! i)e"'(^-2)V4^^^(_l) ^ exponentially smaU terms . (B.91) 

32 



Collecting together (B.89) and (B.91), we conclude that ^^^2(7) exponentially 
small as 7 J, 0. 

Next let us discuss the case N = 2. Here C^. 211+2 — 2, 211-2 — ^ ^ ^ 2, 
C2; 2,1,1 = for / 7^ 1, and C2; 2,1,1 = 2. Thus 



^2:2(7) 



1 1 



1. 



2 - 4^""' - 7e - 2^2,^;2(7) 



(B.92) 



so that, using (B.81), we get that ^2 2(7) exponentially small. 
Finally, for iV = 4, using (2.41), we can write 



E(-i)' 



/=0 



^4; 2,1. 



1+2 



2,1,1-2 



-7A4 



while, from (2.42), we ffet 

00 



l^-' + ^^4,.;l(7) 



-3 + 3e-"'Z4-(7) . 



(B.93) 



(B.94) 



As Z4 (7) is exponentially small, it follows from (B.79) that 2(7) exponentially 
small. 
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C Large- TV Limit 



In this Appendix we want to discuss the iV — )• oo hmit of the finite-size-scahng 
functions for the one-dimensional iV-vector universahty class. We will first discuss 
the derivation using the standard large- iV formalism; then we will show, in two 
particular cases (the spin-1 and spin-2 susceptibilities), how to recover these results 
through a direct evaluation of the iV — )• oo limit of the expressions given in Section 
4.2.2 . 



C.l Review of Results from Standard Large- Formalism 

Let us thus start with the standard large- iV formalism [45]. Consider, on a one- 
dimensional lattice of length L with periodic boundary conditions, the standard 
iV-vector Hamiltonian 

n{{(T}) = -jJ2a,-a,+^ (C.l) 

X 

and the partition function 

Z = J Do-e-^"*^" . (C.2) 

As is well known, the iV — )• oo limit must be taken with J/N fixed. We will therefore 
introduce a rescaled coupling J = J/N. It then turns out [45] that in the iV — )• oo 
limit all correlation functions can be expressed in terms of a mass parameter niL 
which is related to the coupling J by the gap equation 

where p = 27rn/L, the sum runs over n = 0, . . . — 1, and p = 2sin(p/2). The 
summation in (C.3) can be performed exactly, and thus one gets 

J = , coth (^L arcsinh — . (C.4) 



mL\j4: + ml ^ 2 / 

We can now take the limit iV — )• oo at fixed J = J/N and fixed L. All the two- 
point correlation functions (and indeed all the correlation functions) can be easily 
computed [45]: the result is 

G„M^.J:L) = (jjE-^j (C.5a) 

^ cosh" 1(1 -2.T)arcsmhm,/2] fo,, „ < ^. ^ ^ (c.6b) 

cosh [L arcsinhm/,/2] 
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From this expression one immediately gets for the susceptibihties: 



Xoo,k{J ; L) 



1 



1 



cosh*^ [{L — 2a;) arcsinh 

x=0 



cosh*^ [L arcsinh rriL / 2] 
For A; = 1,2 one gets simpler expressions: 

1 



Xoo,i{J]L) 



I d_ 



^—cothfL arcsinh^) 

mL\/4: + ml ^ 2 / 



(C.6a) 
(C.6b) 

(C.7) 
(C.8) 



Analogously one can compute the correlation lengths. For example, in the spin-1 
channel we get 

Ci\j-^L) = — . (C.9) 

rriL 

Having taken the limit N — )• oo, J — )• oo at fixed J = J/N and fixed L, we can 
now take either one of two further limits: 

(a) The standard infinite -volume limit L — )• oo at fixed J. This limit is trivial and 
corresponds simply to the substitution of all sums by the corresponding integrals 
and the parameter niL by moo- In particular the gap equation becomes 

dp 1 1 



J 



TT 27r p2 + m 



4 + m 



(C.IO) 



(b) The finite -size -scaling limit L — )• oo, J — )• oo [hence — )• ooj at fixed S,/L. 
From (C.9) we immediately see that this corresponds to considering the limit L — > 
oo, niL — ^ with itilL = p fixed. The variable p is the natural one in the finite- 
size-scaling limit of the iV = oo model, and all the finite-size-scaling functions will 
be expressed in terms of it. 

Let us first derive the FSS function for the correlation length i^^"/^ Equating 
the right-hand sides of (C.4) and (C.IO) and taking the limit moo — ^ 0, niL — ^ 0, 
L — )• oo with p fixed, we get 



rrir 



tanh 



P 



(C.ll) 



All the other FSS functions can be computed analogously. For the spin-1 and spin-2 
susceptibilities we get 

Xoo,i(L) 



2 P 



Xoo,l(oo) 

XooML) 

Xoo,2(oo) 



tanh 



,2P , P sinhp/2 

tann — \- 5 — — 

2 2 coshV/2 



(C.12) 
(C.13) 
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C.2 Alternate Derivation from Hyperspherical-Harmonics 
Formalism 



Let us now compare our results with those of Section 4.2.2. In that section we 
took first the finite-size-scahng hmit L — )• oo, J — )• oo at fixed 7 = LK{ J) L/(2J) 
and fixed N; now we want to take the further limit N — )• 00, 7 — )• at fixed 
7 = N^}^ We want to show that we recover the same results as in the preceding 
subsection; in other words, we want to show that the two limits commute. 

We need first to find the relation between 7 and p. This is easily obtained if one 
considers in (C.4) the limit J — )• 00, L — )• 00, niL — ^ with 7 = L/(2J) and p fixed. 
We get 

4 = - coth^ . (C.14) 
7 p 2 

Let us begin by computing the limit iV — )• 00, 7 — )• at fixed 7 = of the 
partition-function scaling function 

00 

zPil) = J2^N,ie-^'--' (C.15) 
/=o 

[cf. (4.77)]. Since 7 is tending to zero, it is natural to apply generalized Poisson 
summation formulae of Appendix B. Using (B.28c) we get 

= -5i.(f)"'^^e.p[,(,Y-2)V(4A0]x 

J (/tte-^*'/(4"') f-2zsin|j (C.16) 

— oo-t-«'e 

where e > is arbitrary (the integral is independent of e). The large- iV asymptotic 
expansion of this integral can be obtained by the standard saddle-point technique. 
We rewrite the integral as 

+ oo+ie 2 

t/tt(^-2zsin0 e^^W (C.17) 

— oo+«'e 

with 

m = log(-2zsin0 . (CIS) 

We must now find a saddle point, i.e. a solution of f'{t) = with Imt > 0. We 
find t = ip where p is the unique positive solution of (C.14). Expanding around the 
saddle point we get 

+ oo+ie 2 

dtt (-2zsin0 e^^W = 

-oo+«'e 

1/2 

z,(2sinh0%^/(-) (-f [l + 0(l/iV)] . (C.19) 



^This is clearly the correct scaling, since 7 L/(2J) = (L/2NJ). 
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Collecting everything together we get 



2 ' 



1 



exp 



47 ' 8sinhV/2 

47 



-1/2 



+ 



2 sinh — 



X 

2-JV 



[1 + 0(l/iV)] (C.20) 



In order to compute the large- iV behavior of the finite-size-scaling functions for 
the susceptibilities, we must also evaluate the large- iV behavior of more general 
sums of the form 

00 



(C.21) 



k=0 



where i? is a rational function of the form 



R{x) 



P{x) 



(C.22) 



and P{x) is a polynomial. The coefficients of P{x) and coefficients fii are in general 
iV-dependent. These series can be handled using the generalized Poisson summation 
formula (B.50). As an example let us determine the large- iV behavior of the sum 



00 



_^[fc+(jV_2)/2)]ViV 



(C.23) 



The iV-dependence of the denominator is the one which appears in the finite-size- 
scaling functions of the susceptibilities. We want to compute its large- iV behavior 
for a and 7 fixed. Using (B.50) with f{x) = exp[— (7/iV)(a; — {N — 2)/2)^], we can 
rewrite the sum as 



1 UN 
27r I 7 



1/2 



+ oo+«e 



oo+«e 



^^g.(iV-2)t/2g-iVtV(47) X 



,28* . 



N 



/Ja + 2 + — ;e 



(C.24) 

where fq(l3;z) is defined in (B.54). In order to compute the limit iV — )• 00 of the 
integral we must discuss the large- iV expansion of fq(a-\-N/2; z). The leading order 
is easily obtained if one notices that it is independent of a. In this case one can use 
the recursion relation (B.55) to compute it. For q = 1 we get for iV — )• 00 



/,(« + f;.) = |(1 



+ " + f;*) + 0(11 N'') . 



from which 



fi{a + —-,2 



2{l-z) 



l-N 



[1 + 0(1/N)] . 



N 1 + z 

A similar formula, which can be proved by induction, is valid for generic q: 



/g a + 



N 



-; z 



N 



(1 + ^)'^ 



q-N 



[1 + 0(1/N)] . 



(C.25) 
(C.26) 

(C.27) 
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Using this expansion one can rewrite (C.24) in the hmit iV — )• oo as 

1/2 



I ttN 



2 ^ " 



+ oo+«e 



dt sint 



1 



oo+«e 



1 + 



(C.28) 



where /(t) is given in (CIS). The large- iV expansion of the remaining integral is 
then obtained using the same method used for the partition function. We get finally 



1 



N,k 



k=0 



(k + N/2 + a)i 



-7[fc+(iV-2)/2]ViV 



(C.29) 



7 



-1/2 



N 



exp 



+ 



47 8sinh>/2 



-1/2 



sinh p tanh ■ 



2 sinh ■ 



-N 



[1 + 0(1/N)] . 



(C.30) 



Generic sums of the type (C.21) can be handled exactly in the same way using 
the generalized Poisson formula (B.50). In this case what one needs is the large- 
N behavior of the kernel Fn-r(z). To get explicit formulae we must specify the 
iV-dependence of the coefficients jSi. We will assume jSi = N/2 -\- which is 
the dependence of the sums appearing in the finite-size-scaling functions of the 
susceptibilities. Using the fact that the generic kernel is obtained by summing and 
taking derivatives with respect of z of /g, using (C.27) we see that generically the 
large- iV behavior of Fn-r(z) is given by 



Fn;r{^) = N^il - z)-''Fn{z) [1 + 0(1 /N)] 



(C.31) 



where Fr{z) is a rational function of z independent of iV, and p is an integer. Then 
we obtain the general formula 



Y,J^N,kR{k)e 

k=0 



-f[k+{N-2)/2]yN 



7 



-l/2jyp 



exp 



47 



1 



+ 



1 



-1/2 



sinh pFR(e 2sinh 



47 ' 8sinh>/2_ 

N 

[1 + 0(1/N)] 



P 



(C.32) 



Thus the whole computation boils down to the computation of Fr{z). 

In some cases it is possible to simplify the computation by using a relation 
between the large- iV behavior of different series. Indeed let us differentiate (C.32) 
with respect to 7. Keeping only the leading- iV contributions we get 



1 ^ f N - 

-j^i:^N,kR(k) [k+^ 



_^[fc+(JV_2)/2]ViV 



■ 472 



J2^N,kR{k)e- 



_^[fc+(JV_2)/2]ViV 



[l + 0(l/iV)] , (C.33) 



\k=0 
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where we have used the relation (C.14) to ehminate the terms proportional to dpi d^. 
In the large- iV limit we can of course substitute {k-\-^^Y^Y by {k -\- ^ -\- a){k -\- ^ + 
where a and fi are arbitrary. Thus we get the relation 

£ M^,, R(k) fk + ^ + a)(k + ^ + ^) ,-~.l^HN-2mVN ^ 

k=0 



[1 + 0(l/iV)] 



(C.34) 



This formula will allow us to compute the large- iV behavior of the sums appearing 
in the numerators of the spin-1 and spin-2 susceptibilities. 
Let us start with the spin-1 case. Using (2.39) we get 



N{N - 3) 



l.m 



Using (C.34) with a = -1/2, /3 = -3/2 we get 

-1>-N,l ' - ^ ') 



^ ^^'^ [/ + iV/2- l/2][/ + iV/2-3/2] 



Irri 



-\ 4°^(7) [1 + 0(l/iV)] 



(C.35) 



(C.36) 



and thus 



4t,i(7) 



7 

.P. 



tanh 



2 P 



(C.37) 



which coincides with (C.12). 

In order to evaluate the large- iV limit of the finite-size-scaling function for the 
spin-2 susceptibility we need to evaluate the series (B.67) and (B.68). Using again 
(C.34) we get 



1 



fc=o ((A; + (iV-2)/2)'-l)- 



-l[k+(N-2)l2flN 



'^)'' ff:Ar^.,,e-['=+(^-^)/^]V^') [1 + 0(l/iV)] 



(C.38) 



Vfc=0 



Then using (B.86), (B.87) and (B.90) we get 



l,m=0 

and thus 



A 



JV;/,i 



N 

T 

(0) 



tanh^ — -\-—- 



cos 



hV/2 



4°)(7) [1 + 0(l/iV)] (C.39) 



tanh^^ + ^- ' 



V;oo,2V/y — 2 ' 2coshV/2 
Using (C.14), one immediately sees that (C.40) agrees with (C.13). 



(C.40) 
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